AFIT/GAE/ENY /96M-03 


COMPRESSIBLE  TURBULENCE  MEASUREMENTS  IN  A  SUPERSONIC 
BOUNDARY  LAYER  WITH  IMPINGING  SHOCK  WAVE  INTERACTION 


THESIS 

Michael  J.  Meyer 
First  Lieutenant,  USAF 

AFIT/GAE/ENY/96M-03 


Modal  deformation 


MAF02601 


March  1996 


Master’s  Thesis 


COMPRESSIBLE  TURBULENCE  MEASUREMENTS  IN  A  SUPER¬ 
SONIC  BOUNDARY  LAYER  WITH  IMPINGING  SHOCK  WAVE  INTER¬ 
ACTION 

Michael  J.  Meyer 


Air  Force  Institute  of  Technology,  WPAFB  OH  45433-6583 


AFIT/GAE/ENY/96M-03 


Dr.  James  McMichael 
AFOSR/NA,  BOLLING  AFB,  MD 


Distribution  Unlimited 


This  study  used  laser  Doppler  velocimetry  techniques  to  measure  the  turbulence  intensities  and  Reynolds  shear 
stresses  in  a  Mach  2.9  turbulent  flat  plate  with  an  impinging  shock  wave.  The  shock  strength  was  varied, 
using  three  wedge  shock  generators  with  turning  angles  of  5°,  7°,  and  10".  Each  flow  field  was  studied  at  four 
locations,  X  =  — 1.32<50, 06o,  1.32 60, 2.64 6a,  relative  to  the  shock  wave/boundary  layer  interaction.  An  additional 
streamwise  profile,  parallel  to  the  wall  at  a  height  of  j-  =  0.66.  This  study  showed  that  the  velocity  fluctuations 

)  were  energized  by  the  shock  wave/boundary  layer  interaction,  and  while  remaining  at  a  higher  value 
than  the  upstream  stations,  decreased  in  value  with  distance  downstream  of  the  interaction.  The  Reynolds  shear 
stress  — pu'v '  followed  a  similar  pattern,  sharply  increasing  immediately  after  the  interaction  and  dissipating 
with  distance  downstream  of  the  interaction. 


Shock  Waves,  Supersonic,  Wind  Tunnel,  Turbulence,  Laser  Doppler  Velocimetry, 
Shock  Wave/Boundary  Layer  Interaction 


UNCLASSIFIED 


UNCLASSIFIED 


UNCLASSIFIED 


UL 


✓ 


The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official  policy  or 
position  of  the  Department  of  Defense  or  the  U.  S.  Government. 


AFIT/GAE/ENY /96M-03 


COMPRESSIBLE  TURBULENCE  MEASUREMENTS  IN  A  SUPERSONIC 
BOUNDARY  LAYER  WITH  IMPINGING  SHOCK  WAVE  INTERACTION 


THESIS 


Presented  to  the  Faculty  of  the  Graduate  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 
In  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Master  of  Science  in  Aeronautical  Engineering 


Michael  J.  Meyer,  B.S.E.,  Aerospace  Engineering 
First  Lieutenant,  USAF 


March,  1996 


Approved  for  public  release;  distribution  unlimited 


Acknowledgements 


I  am  deeply  indepted  to  my  thesis  advisor,  Major  Thomas  A.  Buter.  His  guidance  and 
knowledge  were  critical  to  this  research.  His  example  and  motivation  greatly  inspired  me  to  finish 
this  project.  I  would  not  be  here  if  it  was  not  for  him.  Also,  Dr  Rodney  Bowersox  has  provided 
invaluable  expertise  and  encouragement  throughout  this  project.  In  addition,  I  am  grateful  to  Lt 
Joel  Luker  and  Lt  Chad  Hale  for  allowing  me  to  learn  the  laboratory  systems  while  they  were 
trying  to  finish  their  own  thesis.  I  will  not  forget  their  help  and  friendship.  I  am  also  thankful  to 
my  supervisor  at  Wright  Laboratory,  Mr.  Edward  O.  Roberts.  He  generously  allowed  me  the  time 
to  pursue  this  project  and  always  encouraged  me  throughout  my  Master’s  Degree.  Finally,  I  would 
like  to  thank  my  wife  Laura  for  her  patience  and  understanding.  I  have  lost  count  of  the  times  she 
has  rebuilt  my  determination  through  her  love  and  encouragement. 


Michael  J.  Meyer 


Table  of  Contents 


Page 

Acknowledgements .  ii 

List  of  Figures  .  vi 

List  of  Tables .  xii 

Abstract .  xiii 

I.  Introduction .  1-1 

1.1  Introduction .  1-1 

1.2  Thesis  Objective .  1-4 

1.3  Document  Roadmap .  1-4 

II.  Governing  Equations .  2-1 

2.1  Forms  of  the  Governing  Equations .  2-1 

2.2  Compressible  Navier-Stokes  Equations .  2-1 

2.3  Reynolds-Averaged  Navier-Stokes  Equations .  2-3 

2.4  Favre-Averaged  Navier-Stokes  Equations  .  2-5 

III.  Experimental  Apparatus  and  Procedures .  3-1 

3.1  Facility  Description  .  3-1 

3.2  Description  of  Wind  Tunnel .  3-2 

3.3  Laser  Doppler  Velocimeter  .  3-4 

3.4  Photographic  Analysis  Apparatus .  3-9 

3.5  Pressure  and  Temperature  Measurement  .  3-9 

3.6  Model  Description .  3-9 

3.7  Wind  Tunnel  Test  Procedure .  3-11 

3.8  Traverse  and  Laser  Alignment  .  3-12 

3.9  Burst  Spectrum  Analyzer  Settings .  3-13 

iii 


Page 

IV.  Data  Reduction  Techniques .  4-1 

4.1  Velocity  Calculation .  4-1 

4.2  Separation  of  Mean  and  Turbulent  Mach  Numbers .  4-1 

4.3  Mach  Number  Calculation  .  4-2 

4.4  Density  and  Temperature  Calculation .  4-2 

4.4.1  Flat  Plate .  4-3 

4.5  Boundary  Layer,  Displacement  and  Momentum  Thickness  Calculations  4-3 

4.6  Turbulence  Statistics .  4-4 

4.7  Wall  Shear  Stress  .  4-5 

4.8  Turbulent  Kinetic  Energy .  4-6 

V.  Results  and  Discussion .  5-1 

5.1  Roadmap .  5-1 

5.2  Comments  on  Coordinate  Systems .  5-1 

5.3  Test  Conditions  .  5-2 

5.4  Flow  Visualization .  5-2 

5.4.1  Flat  Plate .  5-2 

5.5  Effective  Flow  Turning  Angle .  5-4 

5.5.1  Shock  Wave/Boundary  Layer  Interaction  .  5-5 

5.6  Pressure  and  Temperature  Measurements .  5-12 

5.7  Flat  Plate  LDV  Results .  5-14 

5.7.1  Flat  Plate  Mean  Flow  Characteristics .  5-15 

5.7.2  Flat  Plate  Velocity  Fluctuations .  5-17 

5.7.3  ’’Incompressible”  Reynolds  Shear  Stress .  5-20 

5.8  Shock  Wave/Boundary  Layer  Interaction  LDV  Results .  5-21 

5.8.1  Scaling  Factors .  5-21 

5.8.2  Comparison  to  Rose  and  Johnson  Shock  Wave/Boundary  Layer 

Interaction  .  5-21 


IV 


Page 

5.8.3  Shock  Wave  Boundary  Layer  Results .  5-28 

5.9  Effect  of  Shock  Strength .  5-32 

5.9.1  Fluctuating  Mach  Number  .  5-38 

5.9.2  Wall  Shear  Stress .  5-38 

VI.  Conclusions  and  Recommendations .  6-1 

6.1  Summary  and  Conclusions  .  6-1 

6.2  Recommendations .  6-2 

Appendix  A.  Rose  and  Johnson  1975  Data .  A-l 

A.l  Rose  and  Johnson  1975  Tabular  Data .  A-l 

A. 2  Rose  and  Johnson  1975  Graphical  Data .  A-4 

Appendix  B.  Error  Analysis .  B-l 

Appendix  C.  Turbulence  Statistics  .  C-l 

C.l  Skewness .  C-l 

C.2  Flatness .  C-l 

C.3  Intermittency .  C-17 

C. 4  Velocity  Correlation  Coefficient .  C-17 

Appendix  D.  Tabulated  Data .  D-l 

D. l  Flat  Plate  Tabulated  Data .  D-l 

D.2  5°  Wedge  (Oeff  —  8.2°)  Tabulated  Data .  D-2 

D.3  7°  Wedge  (6eff  =  10.4°)  Flow  Tabulated  Data .  D-9 

D.4  10°  Wedge  ( 0e//  =  13.5°)  Flow  Tabulated  Data .  D-17 

Bibliography .  BIB-1 

Vita .  VITA-1 


v 


List  of  Figures 

Figure  Page 

3.1.  Tunnel  Schematic .  3-1 

3.2.  AFIT  Mach  3  Wind  Tunnel .  3-1 

3.3.  Seeder  System  Schematic .  3-2 

3.4.  Photograph  of  Test  Section  in  AFIT  Mach  3  Wind  Tunnel .  3-3 

3.5.  Vacuum  Pumps  for  AFIT  Mach  3  Wind  Tunnel .  3-4 

3.6.  Laser  and  Photomultiplier  Tubes .  3-6 

3.7.  Laser  Optics  and  Traverse .  3-7 

3.8.  Data  Acquisition  Computer  and  Burst  Spectrum  Analyzers  .  3-8 

3.9.  Shadowgraph  Setup .  3-9 

3.10.  Schlieren  Setup .  3-10 

3.11.  Photograph  of  Shock-Generating  Wedge  and  Tunnel  Floor .  3-10 

3.12.  Pressure  Tap  Locations .  3-11 

3.13.  Laser  and  Traverse  Alignment  Relative  to  Test  Section .  3-12 

5.1.  Wind  Tunnel  and  Shock  Wave/Boundary  Layer  Interaction  Coordinate  Systems  5-2 

5.2.  Schlieren  of  Flat  Plate  Using  Horizontal  Knife  Edge .  5-3 

5.3.  Shadowgraph  of  Flat  Plate  Flowfield .  5-4 

5.4.  Schlieren  of  5°  Wedge  Flow  Interaction  Using  Horizontal  Knife  Edge  (Forward 

Window) .  5-6 

5.5.  Shadowgraph  of  5°  Wedge  Flow  Interaction  Flowfield  (Forward  Window)  ....  5-7 

5.6.  Schlieren  of  5°  Wedge  Flow  Interaction  Using  Horizontal  Knife  Edge  (Rear  Win¬ 
dow)  .  5-8 

5.7.  Shadowgraph  of  5°  Wedge  Flow  Interaction  Flowfield  (Rear  Window) .  5-9 

5.8.  Schlieren  of  7°  Wedge  Flow  Interaction  Using  Horizontal  Knife  Edge .  5-10 

5.9.  Shadowgraph  of  7°  Wedge  Flow  Interaction  Flowfield .  5-11 

5.10.  Schlieren  of  10°  Wedge  Flow  Interaction  Using  Horizontal  Knife  Edge .  5-12 

vi 


Figure  Page 

5.11.  Shadowgraph  of  10°  Wedge  Flow  Interaction  Flowfield .  5-13 

5.12.  Pressure  Distribution  over  Flat  Plate  with  Impinging  Shock  Wave .  5-14 

5.13.  Pressure  Distribution  over  Shock  Generating  Wedges .  5-15 

5.14.  Normalized  Flat  Plate  Velocity  Profile .  5-16 

5.15.  Flat  Plate  Mach  Number  Profile .  5-17 

5.16.  Flat  Plate  Fluctuating  Mach  Number  Profile .  5-18 

5.17.  Flat  Plate  Streamwise  Velocity  Fluctuations  Scaled  by  Local  Velocity .  5-18 

5.18.  Flat  Plate  Streamwise  Velocity  Fluctuations  Scaled  by  Edge  Velocity .  5-19 

5.19.  Flat  Plate  Cross-stream  Velocity  Fluctuation  Scaled  by  Local  Velocity .  5-19 

5.20.  Flat  Plate  Cross-stream  Velocity  Fluctuation  Scaled  by  Edge  Velocity  .  5-20 

5.21.  Incompressible  Reynolds  Shear  Stress  Scaled  by  Local  Velocity .  5-21 

5.22.  Normalized  Velocity  Profile  Comparison  at  Upstream  Stations .  5-22 

5.23.  Normalized  Velocity  Profile  Comparison  at  First  Downstream  Stations .  5-23 

5.24.  Normalized  Velocity  Profile  Comparison  at  Second  Downstream  Stations  ....  5-23 

5.25.  Streamwise  Velocity  Fluctuation  Comparison  at  the  Upstream  Station .  5-24 

5.26.  Streamwise  Velocity  Fluctuation  Comparison  at  \st  Downstream  Stations  .  .  .  5-25 

5.27.  Streamwise  Velocity  Fluctuation  Comparison  at  the  2nd  Downstream  Stations  .  5-25 

5.28.  Cross-stream  Velocity  Fluctuations  Comparison  at  the  Upstream  Station  ....  5-26 

5.29.  Cross-stream  Velocity  Fluctuation  Comparison  at  1'*  Downstream  Station  .  .  .  5-27 

5.30.  Cross-stream  Velocity  Fluctuation  Comparison  at  2nd  Downstream  Station  .  .  .  5-27 

5.31.  Reynolds  Shear  Stress  Scaled  by  Wall  Shear  Stress  for  the  Upstream  Station  .  .  5-28 

5.32.  Reynolds  Shear  Stress  Scaled  by  Wall  Shear  Stress  for  2nd  Downstream  Station  5-29 

5.33.  5°  Wedge  (&eff  =  8.2°)  Velocity  Profiles  at  Each  Measurement  Station .  5-30 

5.34.  5°  Wedge  (0e//  =  8.2°)  Mach  Number  Profiles  at  each  Measurement  Station  .  .  5-30 

5.35.  5°  Wedge  (0eff  =  8.2°)  Streamwise  Velocity  Fluctuations .  5-31 

5.36.  5°  Wedge  (6efj  =  8.2°)  Cross-stream  Velocity  Fluctuations .  5-32 

5.37.  5°  Wedge  (Oeff  =  8.2°)  Incompressible  Reynolds  Shear  Stress  Scaled  by  Local 

Velocity  .  5-33 

vii 


Figure  Page 

5.38.  5°  Wedge  (8eff  =  8.2°)  Turbulent  Kinetic  Energy .  5-33 

5.39.  Streamwise  Velocity  Profile  at  Height  from  Wall  V  =  0.666o .  5-34 

5.40.  Streamwise  Mach  Number  Profiles  at  Height  from  Wall  Y  —  0.66 80  .  5-35 

5.41.  Streamwise  Velocity  Fluctuations  at  Height  from  Wall  Y  —  0.666,, .  5-36 

5.42.  Cross-stream  Velocity  Fluctuations  at  Height  from  Wall  Y  =  0.666o .  5-36 

5.43.  Incompressible  Reynolds  Shear  Stess  Scaled  by  Local  Velocity  at  Height  from  Wall 

Y  =  0.666o .  5-37 

5.44.  Turbulent  Kinetic  Energy  at  Height  from  Wall  Y  =  0.666o .  5-38 

5.45.  Comparisons  of  Wall  Shear  Stress .  5-40 

5.46.  Skin  Friction  Coefficients  for  the  Shock  Wave/Boundary  Layer  Interaction  .  .  .  5-41 

A.l.  Mean  Velocity  Profiles .  A-4 

A. 2.  Streamwise  Velocity  Fluctuations .  A-4 

A.3.  Crosswise  Velocity  Fluctuations .  A-5 

A.4.  Reynolds  Normal  Stress  .  A-5 

A.5.  Reynolds  Shear  Stress  .  A-6 

C.l.  X  =  —1.32 S0  Streamwise  Skewness .  C-l 

C.2.  X  =  — 1.3260  Cross-stream  Skewness .  C-2 

C.3.  X  =  06o  Streamwise  Skewness .  C-2 

C.4.  X  =  06o  Cross-stream  Skewness .  C-3 

C.5.  X  =  1.3260  Streamwise  Skewness .  C-3 

C.6.  X  =  1.3260  Cross-stream  Skewness .  C-4 

C.7.  X  =  2.6460  Streamwise  Skewness .  C-4 

C.8.  X  =  2.64 S0  Cross-stream  Skewness .  C-5 

C.9.  Streamwise  Skewness  at  a  Height  from  Wall  Y  =  0.666o .  C-5 

C.10.  Cross-stream  Skewness  at  a  Height  from  Wall  Y  =  O.6660 .  C-6 

C.ll.  5°  Wedge  Streamwise  Skewness  at  each  Station .  C-6 

C.12.  5°  Wedge  Cross-stream  Skewness  at  each  Station .  C-7 

viii 


Figure  Page 

C.13.  7°  Wedge  Streamwise  Skewness  at  each  Station .  C-7 

C.14.  7°  Wedge  Cross-stream  Skewness  at  each  Station .  C-8 

C.15.  10°  Wedge  Streamwise  Skewness  at  each  Station .  C-8 

C.16.  10°  Wedge  Cross-stream  Skewness  at  each  Station .  C-9 

C.17.  X  =  — 1.3260  Streamwise  Flatness .  C-9 

C.18.  X  =  —1.32 60  Cross-stream  Flatness .  C-10 

C.19.  X  =  0 S0  Streamwise  Flatness  .  C-10 

C.20.  X  =  0 S0  Cross-stream  Flatness  .  C-ll 

C.21.  X  =  1.32 S0  Streamwise  Flatness .  C-ll 

C.22.  X  =  1.32 60  Cross-stream  Flatness .  C-12 

C.23.  X  =  2.64 60  Streamwise  Flatness .  C-12 

C.24.  X  =  2.64 80  Cross-stream  Flatness .  C-13 

C.25.  Streamwise  Flatness  at  a  Height  from  Wall  Y  =  O.66<50 .  C-13 

C.26.  Cross-stream  Flatness  at  a  Height  from  Wall  Y  —  0.66<5o .  C-14 

C.27.  5°  Wedge  Streamwise  Flatness  at  each  Station .  C-14 

C.28.  5°  Wedge  Cross-stream  Skewness  at  each  Station .  C-15 

C.29.  7°  Wedge  Streamwise  Flatness  at  each  Station .  C-15 

C.30.  7°  Wedge  Cross-stream  Flatness  at  each  Station .  C-16 

C.31.  10°  Wedge  Streamwise  Flatness  at  each  Station .  C-16 

C.32.  10°  Wedge  Cross-stream  Flatness  at  each  Station .  C-17 

C.33.  X  =  — 1.32<50  Intermittency  Function .  C-18 

C.34.  X  =  06o  Intermittency  Function .  C-18 

C.35.  X  =  1.32<50  Intermittency  Function .  C-19 

C.36.  X  =  2.6460  Intermittency  Function .  C-19 

C.37.  5°  Wedge  Intermittency  Function .  C-20 

C.38.  7°  Wedge  Intermittency  Function .  C-20 

C.39.  10°  Wedge  Intermittency  Function  .  C-21 


IX 


Figure  Page 

C.40.  Velocity  Correlation  Coefficient  at  X  =  —1.32 60  .  C-21 

C.41.  Velocity  Correlation  Coefficient  at  X  =  0 60  .  C-22 

C.42.  Velocity  Correlation  Coefficient  at  X  =  1.32£e .  C-22 

C.43.  Velocity  Correlation  Coefficient  at  X  =  2.6460 .  C-23 

C.44.  5°  Wedge  Correlation  Coefficient  .  C-23 

C.45.  7°  Wedge  Correlation  Coefficient  .  C-24 

C. 46.  10°  Wedge  Correlation  Coefficient .  C-24 

D. l.  Tabulated  Data  for  Flat  Plate  at  x  =  48  cm .  D-l 

D.2.  5°  Wedge  (8eff  =  8.2°)  Flow  X  =  —1.32 80  Coincident  Data .  D-2 

D.3.  5°  Wedge  (8eff  =  8.2°)  Flow  X  —  0 60  Coincident  Data .  D-3 

D.4.  5°  Wedge  ( 6ejj  —  8.2°)  Flow  X  —  1.3260  Coincident  Data .  D-4 

D.5.  5°  Wedge  {Qe}j  —  8.2°)  Flow  X  =  2. 6460  Coincident  Data .  D-5 

D.6.  5°  Wedge  (de//  =  8.2°)  Flow  X  =  2.64 60  Coincident  Data  (Cont.) .  D-6 

D.7.  5°  Wedge  ( 0e//  =  8.2°)  Flow  Y  =  0.665o  Coincident  Data .  D-7 

D.8.  5°  Wedge  ( 0e//  =  8.2°)  Flow  Y  =  0.666o  Coincident  Data  (Cont.) .  D-8 

D.9.  7°  Wedge  (de//  =  10.4°)  Flow  X  =  —1.32 60  Coincident  Data .  D-9 

D.10.  7°  Wedge  (Oefj  =  10.4°)  Flow  X  =  0 60  Coincident  Data .  D-10 

D.ll.  7°  Wedge  (8eff  =  10.4°)  Flow  X  =  1.32 S0  Coincident  Data .  D-ll 

D.12.  7°  Wedge  (8eff  =  10.4°)  Flow  X  =  1.3260  Coincident  Data  (Cont.) .  D-12 

D.13.  7°  Wedge  (0e//  =  10.4°)  Flow  X  —  2.64<50  Coincident  Data .  D-13 

D.14.  7°  Wedge  ( 0e//  =  10.4°)  Flow  X  —  2.64§0  Coincident  Data  (Cont.) .  D-14 

D.15.  7°  Wedge  (9eff  =  10.4°)  Flow  Y  =  0.666o  Coincident  Data .  D-15 

D.16.  7°  Wedge  (8ejf  =  10.4°)  Flow  Y  —  0.666o  Coincident  Data  (Cont.) .  D-16 

D.17.  10°  Wedge  (8eff  —  13.5°)  Flow  X  =  —1.32 60  Coincident  Data .  D-17 

D.18.  10°  Wedge  ( 8e/f  —  13.5°)  Flow  X  —  0 60  Coincident  Data .  D-18 

D.19.  10°  Wedge  {8efs  =  13.5°)  Flow  X  =  1.32 80  Coincident  Data .  D-19 

D.20.  10°  Wedge  (8eff  =  13.5°)  Flow  X  =  1.32 60  Coincident  Data  (Cont.)  .  D-20 


x 


Figure  Page 

D.21.  10°  Wedge  (0eff  =  13.5°)  Flow  X  =  2.64 S0  Coincident  Data .  D-21 

D.22.  10°  Wedge  (0ejj  —  13.5°)  Flow  X  —  2.64 60  Coincident  Data  (Cont.)  .  D-22 

D.23.  10°  Wedge  {0efS  =  13.5°)  Flow  Y  =  0.66do  Coincident  Data .  D-23 

D.24.  10°  Wedge  (0e//  =  13.5°)  Flow  Y  =  0.665o  Coincident  Data  (Cont.) .  D-24 


xi 


List  of  Tables 

Table  Page 

3.1.  Quick  Menu  Constant  Settings .  3-13 

3.2.  Soft  Menu  Constant  Settings .  3-14 

3.3.  BSA  Program  Menu  Constant  Settings .  3-14 

5.1.  Normalized  Shock  Wave  Interaction  Pressure  Increase .  5-14 

5.2.  Scaling  Factors  for  Present  Study  and  Rose  and  Johnson .  5-21 

5.3.  Downstream  Measurement  Stations  in  Boundary  Layer  Thicknesses  from  Interac¬ 
tion  .  5-22 

5.4.  Pressure  Gradients  in  Pa/m  Used  in  Calculation  of  Wall  Shear  Stress .  5-39 

A.l.  Mean  Velocity  Profiles .  A-l 

A.2.  Streamwise  Velocity  Fluctuations .  A-l 

A.3.  Crosswise  Velocity  Fluctuations .  A-2 

A.4.  Reynold’s  Normal  Stress .  A-2 

A.  5.  Reynold’s  Shear  Stress .  A-3 

B. l.  Error  Summary:  Values  marked  by  *  represent  changes  from  Luker .  B-2 


xii 


AFIT/GAE/ENY/96M-03 


Abstract 

This  study  used  laser  Doppler  velocimetry  techniques  to  measure  the  turbulence  intensities 
and  Reynolds  shear  stresses  in  a  Mach  2.9  turbulent  flat  plate  with  an  impinging  shock  wave.  The 
shock  strength  was  varied,  using  three  wedge  shock  generators  with  turning  angles  of  5°,  7°,  and 
10".  Each  flow  field  was  studied  at  four  locations,  X  =  —1.32 S0, 0 60, 1.32 60, 2.64<50,  relative  to  the 
shock  wave/boundary  layer  interaction.  An  additional  streamwise  profile,  parallel  to  the  wall  at  a 
height  of  —  0.66.  This  study  showed  that  the  velocity  fluctuations  ( were  energized 
by  the  shock  wave/boundary  layer  interaction,  and  while  remaining  at  a  higher  value  than  the 
upstream  stations,  decreased  in  value  with  distance  downstream  of  the  interaction.  The  Reynolds 
shear  stress  -pu'v1  followed  a  similar  pattern,  sharply  increasing  immediately  after  the  interaction 
and  dissipating  with  distance  downstream  of  the  interaction. 


COMPRESSIBLE  TURBULENCE  MEASUREMENTS  IN  A  SUPERSONIC 
BOUNDARY  LAYER  WITH  IMPINGING  SHOCK  WAVE  INTERACTION 


I.  Introduction 

1.1  Introduction 

The  interaction  between  supersonic  boundary  layers  and  impinging  shock  waves  represents 
a  significant  problem  for  the  designers  of  high-speed  air  vehicles.  Impinging  shock  waves  appear 
frequently  in  engine  inlets  and  compression  systems.  The  abrupt  pressure  rise  caused  by  the 
impinging  shock  wave  complicates  the  vehicle  structural  design.  If  the  shock  is  strong  enough,  the 
boundary  layer  may  separate,  causing  higher  heat  transfer  rates  and  burdening  the  designer  with 
the  additional  problem  of  thermal  protection  (19). 

This  problem  appeared  with  a  vengeance  on  October  3,  1967  when  an  X-15  high-speed  re¬ 
search  aircraft,  piloted  by  Pete  Knight,  reached  a  Mach  Number  of  6.72  at  an  altitude  just  over 
100,000  feet.  The  aircraft  carried  a  dummy  ramjet  connected  to  the  lower  surface  of  the  X-15  by 
a  pylon.  During  the  flight,  a  shock  wave  from  the  nacelle  of  the  ramjet  impinged  on  the  pylon, 
burning  a  hole  through  the  surface.  In  addition,  the  pylon  itself  created  a  bow  shock  that  im¬ 
pinged  on  the  lower  surface  of  the  X-15,  causing  local  heating  damage.  The  dummy  ramjet  was 
completely  burned  off  the  pylon  and,  during  separation  from  the  aircraft,  punched  a  hole  in  the 
X-15  that  allowed  heated  external  air  into  the  internal  aircraft  structure.  Knight  was  able  to  land 
the  damaged  X-15  but  this  was  the  worst  case  of  aerodynamic  heating  damage  experienced  during 
the  X-15  Program  (10). 

Flight  testing  or  wind  tunnel  testing  of  hypersonic  and  supersonic  air  vehicle  configurations 
to  identify  these  impinging-shock-wave  problems  can  be  cost-prohibitive  and  possibly  hazardous. 
Using  Computational  Fluid  Dynamics  (CFD),  the  designer  will  have  ample  opportunity  to  test  and 


1-1 


refine  the  design  on  the  ground  at  a  much  lower  cost.  However,  before  this  is  possible,  advanced 
CFD  techniques  must  be  developed  to  handle  the  different  flow  structures  expected  in  high-speed 
flight  (19). 

Directly  solving  the  Navier-Stokes  Equations  for  high  Reynolds  number  flow  is  impractical. 
Luker  (13)  shows  that  for  the  present  test  section  (6.35cm  by  6.35  cm  by  80cm,  M  =  2.9,  Res 
roughly  1.6  X  106,  boundary  layer  thickness  delta  roughly  =  0.009  m)  the  Kolmogoroff  turbulent 
length  scale  is  roughly  2  X  10-7  m.  To  resolve  this  flow  we  would  need  4  X  1017  grid  points,  4  X  1028 
calculations,  and,  at  an  optimistic  one  nanosecond  per  calculation,  1.2  X  1012  years  to  complete 
the  numerical  solution  to  the  Navier-Stokes  solutions.  The  most  advanced  computers  of  the  day 
can  not  make  this  calculation  and  most  would  not  wait  for  it.  This  is  a  relatively  simple  flowfield. 
Direct  calculation  of  low  Reynolds  number  flows  have  been  completed  but  for  most  practical  flows 
aerodynamicists  must  resort  to  approximate  methods.  These  approximate  methods  usually  involve 
averaging  the  Navier-Stokes  equations  and  writing  the  dependent  variables  in  terms  of  a  mean  and 
fluctuation  about  that  mean.  However,  when  the  non-linear  Navier-Stokes  equations  are  cast  in 
this  form,  cross  correlation  fluctuation  terms  appear,  leading  to  fewer  equations  than  unknowns. 
These  terms  cause  the  problem  of  turbulent  closure.  To  achieve  closure  is  the  role  of  the  turbulence 
model. 

The  development  of  these  advanced  CFD  techniques  and  turbulence  models  is  greatly  aided 
by  careful  experiments  performed  to  provide  validation  cases  and  verifying  the  accuracy  of  the 
computational  techniques.  These  experiments  not  only  increase  the  basic  understanding  of  the 
flow  structure  and  fluid  motions  but  can  also  provide  input  for  modeling  and  precise  checks  on 
computational  output.  For  research  codes  and  initial  trials  of  turbulence  models,  building  block 
experiments,  providing  phenomenological  modeling  information,  are  needed  (14). 

To  achieve  closure,  simplified  turbulence  models  are  used.  However,  these  models  are  not  rig¬ 
orously  correct  and  are  not  adequate  for  general  predictions  of  flowfields  around  high-speed  vehicles. 
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Compressible,  supersonic  flows  have  drastically  different  properties  from  the  simple,  low-speed  flows 
initially  studied  and  used  to  build  turbulence  models.  Yet,  many  incompressible  turbulence  models 
are  applied  to  compressible,  high-speed  flow  with  little  or  no  theoretical  justification.  This  method 
has  been  reasonably  successful  for  flat  plate  flow  but  turbulent  quantities  change  greatly  for  com¬ 
plex  geometries  and  large  pressure  gradients(18).  The  complexity  of  turbulence  modeling  demands 
great  expertise  and  large  quantities  of  intellectual  energy  and  must  be  guided  by  careful  experi¬ 
ments.  Upon  the  introduction  of  a  new  turbulence  model,  it  must  be  tested  and  validated  against 
quality  data.  To  adequately  test  these  turbulence  models,  a  broad  range  of  flow  field  experiments 
with  careful  measurement  of  the  turbulent  quantities  is  needed  (17). 

This  research  will  attempt  to  meet  the  eight  necessary  criteria  for  a  CFD  validation  case  as 
spelled  out  by  Settles  and  Dodson  (17): 


Baseline  Applicability  :  Each  experiment  must  involve  Mach  3  or  higher  turbu¬ 
lent  flow  with  a  shock  wave-boundary  layer  interaction. 

Simplicity  :  Each  experiment  must  involve  a  simple  ’’building-block”  geometry 
which  can  be  sufficiently  modeled  by  CFD  methods  without  great  difficulty. 

Specific  Applicability  :  Each  experiment  must  provide  a  useful  test  of  a  turbu¬ 
lence  model. 

Well-defined  Experimental  Boundary  Conditions  :  This  condition  requires 
that  all  incoming  boundary  conditions  be  known  including  the  state  of  the  incoming 
boundary  layer.  Also,  a  claim  of  two-dimensional  flow  must  show  experimental  data  to 
verify  the  extent  of  spanwise  flow. 

Well-defined  Experimental  Error  Bounds  :  Each  experiment  must  be  accom¬ 
panied  by  an  analysis  of  the  accuracy  and  repeatability  of  the  data  to  ensure  proper 
code  validation  can  be  conducted. 

Consistency  Criterion  :  No  mutually  inconsistent  results  can  appear  in  the  ex¬ 
periment. 

Adequate  Documentation  of  Data  :  The  data  from  the  experiment  must  be 
tabulated  and  in  a  machine-readable  form. 

Adequate  Spatial  Resolution  of  Data  :  Each  experiment  must  have  adequate 
spatial  resolution  to  resolve  the  key  features  of  the  flow. 
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Using  these  eight  criteria  to  review  of  the  literature,  Settles  and  Dodson  did  not  find  an 
accepted  experiment  for  impinging  shock  wave  -  boundary  layer  interaction.  Secondary  criteria 
include  turbulence  data,  realistic  test  conditions,  nonintrusive  instrumentation,  and  redundant 
measurements  (17). 

The  only  source  found  that  measured  impinging  shock  wave  -  boundary  layer  interaction 
with  nonintrusive  instrumentation  was  Rose  and  Johnson’s  1975  study.  In  this  study,  Rose  and 
Johnson  only  report  on  the  shock  wave  -  boundary  layer  interaction  for  one  shock  system  caused  by 
a  7  degree  ramp.  However,  the  study  does  compare  Laser  Doppler  Velocimetry  (LDV)  data  with 
hot-wire  anemometry  data  for  this  flow  field  (16). 

1.2  Thesis  Objective 

The  objective  of  this  thesis  is  to  determine  the  effect  of  impinging  shock  wave  with  varied 
strengths  on  the  turbulence  characteristics  of  a  boundary  layer  in  Mach  2.9  flow.  A  secondary 
objective  of  this  thesis  is  to  provide  quality  validation  data  for  turbulence  model  developers  that 
meets  the  eight  criteria  outlined  by  Settles  and  Dodson. 

1.3  Document  Roadmap 

As  a  basis  for  the  determination  of  the  necessary  validation  measurements  for  turbulence 
modeling,  the  approximate  form  of  the  governing  equations  of  compressible,  turbulent  flow  are 
presented  in  Chapter  Two,  showing  the  importance  of  the  velocity  cross-correlation  term,  u'v'.  An 
explanation  of  the  experimental  apparatus  and  test  methodology  used  are  presented  in  Chapter 
Three.  The  data  reduction  techniques  used  to  determine  the  turbulent  properties  of  shock  wave¬ 
boundary  layer  interactions  are  addressed  in  Chapter  Four.  The  results  of  this  study  are  presented 
and  discussed  in  Chapter  Five  culminating  in  a  series  of  conclusions  provided  in  Chapter  Six,  along 
with  a  series  of  recommendations  for  further  research. 
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II.  Governing  Equations 


2.1  Forms  of  the  Governing  Equations 

Computational  study  of  turbulent  flow  is  primarily  conducted  by  numerically  solving  the  time- 
averaged  Navier-Stokes  equations.  This  system  of  equations  is  derived  by  breaking  the  dependent 
variables  into  mean  and  fluctuating  components  and  time-averaging  the  entire  equation  set.  This 
process  leads  to  new  terms  which  appear  as  stress  gradients  and  heat-flux  quantities  which  involve 
products  of  the  fluctuating  quantities.  To  accurately  account  for  these  terms,  the  system  must 
be  closed  with  functions  of  mean  flow  variables.  This  is  done  through  a  turbulence  model.  Two 
types  of  averaging  are  used  and  will  be  discussed,  the  classical  Reynolds  averaging  and  Favre  mass- 
weighted  averaging.  The  two  forms  become  identical  when  the  density  fluctuations  are  neglected 

(4). 

2.2  Compressible  Navier-Stokes  Equations 

The  governing  equations  for  continuum  fluid  flow  are  the  Navier-Stokes  equations  presented 
here  in  compressible  form(4). 

Continuity: 

^  +  p(V.V)  =  0  (2.1) 


Momentum: 


P 


DV 

Dt 


=  pf  +  V  ■  7Tjj 


(2.2) 
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Energy: 


DP  dQ 

m  +~m 


-  V  •  q  +  $ 


where  f  represents  the  body  forces  acting  on  the  fluid  and 


(2.3) 


TTij  =  -P<5ij  +  T\j 


(2.4) 


Tij  =  [I 


duk 

dxk 


(2.5) 


$  = 


dui 
Tij  dxj 


(2.6) 


In  turbulent  flow,  rapid  fluctuations  in  the  important  flow  variables  leads  to  very  small  time 
scales  associated  with  the  flowfield.  Averaging  schemes,  such  as  the  Reynolds  or  Favre  averag¬ 
ing,  were  developed  to  lower  the  computational  work  involved  in  solving  flows  with  small  time 
scales.  This  way  only  large  mean  flow  variations  (compared  to  the  turbulence  time  scales)  are 
captured  using  the  time-dependent  Navier-Stokes  equations.  This  leads  to  a  dramatic  decrease  in 
the  computational  burden. 

In  both  averaging  schemes,  the  flow  variables  are  broken  into  a  mean  component  and  a 
fluctuating  component.  The  instantaneous  variable,  h,  consists  of  the  time  averaged  component, 
h,  and  the  fluctuation  from  that  value,  h' . 


h(x>  V’  z>  0  =  h(x,  V,  z)  +  h'(x,  y,  z,  t ) 


(2.7) 


2-2 


It  should  be  noted  that  in  equation  2.7  h  is  allowed  to  vary  with  time,  thus  the  averaging 
process  admits  time-dependent  behavior.  This  process  requires  that  the  time-averaging  period 
must  be  long  compared  to  the  turbulence  time  scales  but  relatively  short  compared  to  the  overall 
flowfield  time  (4). 

2.3  Reynolds-Averaged  Navier-Stokes  Equations 

When  using  Reynolds  averaging,  any  given  primitive  flow  quantity,  </>,  is  expressed  as  the  sum 
of  a  time-averaged  component,  <j>  and  a  fluctuating  component,  o',  as  follows  (4): 


p  =  p  +  p' 

U{  =  Ui  +  u'i 

p  =  p  +  p' 

T  =  T  +  T' 

E  =  E  +  E' 

H  =  H  +  H' 

P  =  p  -f  p! 

p*  =W  +  p*' 

cv  =  c~v  +  c'v 

cp  =  c;+c'p 

The  time  average  is  defined  by 


where  the  characteristic  time  T  is  long  in  comparison  to  the  cycle  times  of  the  fluctuating  compo¬ 
nent,  but  still  permits  gradual  time-dependent  fluid  motion  for  non-stationary  flows. 


To  obtain  the  mean  flow  equations,  these  expressions  for  the  decomposed  variables  are  sub¬ 
stituted  into  the  governing  equations.  These  instantaneous  equations  are  then  time  averaged, 
expanded,  and  simplified  using  the  rules  given  in  (19).  The  results  are  the  Reynolds-averaged 
Navier-Stokes  equations  given  below  (4) 
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Continuity: 


dp  d(puj  +  p'u'j)  _ 
dt  dxi 


Momentum: 

dpui  +  p'u'j  dpujUj  _  _8p_  d(nj  4-  TIj) 

dt  dxj  dxi  dxj 


Energy: 

d(pe^ +  7K)  djpKuj)  _  djumj+u'^-qj-qj) 

dt  dxj  dxj 


(2.9) 


(2.10) 


(2.11) 


In  addition, 


Tij  =  p 


2  duk' 

3  6ijdxk_ 


2 ,  du'k  ‘ 
3  6ijdxk_ 


(2.12) 


(2.13) 


Tij  =  -pu'iU'j  -  Uip'u'j  -  Ujp'u'i  -  p'u'iU'j 


(2.14) 


= 


,  dT 


(2.15) 


qf  =  -ph'0u'i  +  hop'u'i  +  Uip'h'0  -  p'h'0u'i 


(2.16) 


In  the  above  equations,  ea  and  h0  are  stagnation  reference  quantities. 
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2-4  Favre-Averaged  Navier-Stokes  Equations 


When  employing  Favre  averaging,  any  given  primitive  flow  quantity,  <j>,  is  expressed  as  the 
sum  of  a  mass-weighted-averaged  mean  value  (Favre  averaged),  <j>,  and  a  fluctuating  component, 
<f>" .  The  Favre  averaged  quantity  and  the  Favre  turbulent  fluctuation  of  the  quantity  are  defined 
as  (4): 

4>= ^  and  =  (2.17) 

P 


It  should  be  noted  that  although 


<t>"  ±  0 


p<f>  becomes 


forcing 


p(t>  =  p(<t>  +  <t>") 
=  p<f>  +  p<j>" 
=  p<t>  +  p<f>" 

=  pf+Jr 

=  pcp  +  p<t>" 

p<f>"  =  0 


Also,  the  following  relationships  can  be  derived 


(j,  -  $  =  <f>"  -  <j>'  =  <j,"  =  -ElfL 

P 


The  Favre  averaged  equations,  as  presented  by  Marvin  (14)  (and  corrected  here),  are 


Continuity: 


dp  d(puj) 
dt  dxi 


(2.18) 
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Momentum: 


d(puj)  d(pujUj)  _  _dp_  HTij  ~  Pui  uj  ) 

I  r\  I" 


dx,' 


dxi 


<9x,- 


(2.19) 


Energy: 


dph  i  d(phuj)  dp  t  _ 

I  u j 

Cj  C7E 


<9t  '  <9x,-  dt  '  dxt  U ■*  dx 


„  dp  ,  9(-?j  -  ph"u  •)  _  au; 

- 1 - x - H  Tij  — 1 


dx,- 


dxs 


(2.20) 


These  governing  equations  may  be  supplemented  by  the  following  transport  equations. 


Turbulent  Kinetic  Energy: 


d(pk)  d(piljk)  — Ti—irduk  d(pu^k)  d(u”p)  du ’■  d(u"Tik )  d(u'-) 


dt 


+ 


dx  i 


=  ~Pui  uk 


dxk 


dXj 


4-  v - —  + 

dxi  F  dxi  dxk 


T~ik  ‘ 


dxk 


(2.21) 


Reynolds  stress  transport: 


dt  '  da?,- 


=  «'»;)  w  -  wo  W 


a(pu,  _  a(u'.'p)  _  g(«,-p)  g(«")  ,  g(u")  ,  aCuit'Ti,-)  (2.22) 

dx .-  da:.-  da:*  ‘  *  dx ’■  da:*,  ‘da:.-  ' 


da?i 


.  aO'/ftO  _  9(0  _  9(0 

■”  da:,-  '*•?  dxj  da:,- 


These  equations  have  the  same  form  as  the  incompressible  Reynolds  averaged  equations,  ex¬ 
cept  that  the  Reynolds  stresses  pu”u '■  include  the  density  fluctuations,  which  must  be  accounted 
for  in  some  way.  Compressibility  is  generally  included  by  replacing  the  density  with  the  mean 
density  neglecting  the  terms  which  involve  additional  correlations  arising  from  compressibility.  Re¬ 
searchers  using  the  Favre  forms  of  the  governing  equations  cite  that  avoiding  some  of  the  fluctuating 
density-generated  terms  simplifies  closure  of  the  system.  However  since  much  of  the  experimental 
data  provides  RANS-type  information,  modeling  the  effects  of  compressibility  presents  a  problem. 
The  lack  of  an  explicit  fluctuating  density  component  in  the  turbulent  shear  stresses  is  one  of  the 
specific  shortcomings  of  this  model. 
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The  compressible  FANS  turbulence  terms  are  expressed  as 


mi 


rT. 

« 3 


-pu"u" 


Qi  =  PKui 


(2.23) 


It  can  be  seen  that  the  conservation  and  supplemental  equations  have  a  greater  number  of 
unknowns  than  equations.  This  is  the  problem  of  closure  and  must  be  remedied  through  turbulence 
modeling.  A  detailed  discussion  of  turbulence  modeling  is  beyond  the  scope  of  this  thesis.  The 
reader  is  referred  to  Marvin  for  an  excellent  overview  of  the  direction  of  modeling  efforts  (14). 
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III.  Experimental  Apparatus  and  Procedures 


3.1  Facility  Description 


The  experimental  work  presented  here  was  performed  in  the  AFIT  Supersonic  Wind  Tunnel 


Facility.  A  schematic  of  the  facility  is  shown  in  Figure  3.1  and  Figure  3.2 


High  Pressure  Air 
100  PSI 


Vacuum  System 


Gate  Valve 


Tunnel  Schematic 


- - 

Settling  Chamber 

^Nozzle 

Test  Section 

Test  Section 

Diffuser 

- - — - 

Figure  3.2  AFIT  Mach  3  Wind  Tunnel 
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3.2  Description  of  Wind  Tunnel 


The  wind  tunnel  used  for  this  research  is  a  Mach  2.9  combination  blowdown-drawdown  tunnel, 
shown  in  Figure  3.2  supplied  by  100  psi  high  pressure  shop  air  exhausting  into  a  vacuum  system. 
The  high  pressure  air  is  supplied  by  an  Atlas  Compco  compressor  facility  with  an  approximate  mass 
flow  capacity  of  0.45  kg/sec.  The  system  air  is  then  passed  through  two  Pioneer  Air  Systems  Inc. 
Model  R500A  refrigeration  dryers  and  into  a  centrifugal  particle  and  moisture  separator.  Finally, 
the  air  passes  through  a  fiber-reinforced  paper  filter  before  entering  the  settling  chamber.  The 
compressed  air  is  released  into  the  tunnel  using  an  El-o-Matic  ED10  Ball  Valve.  In  the  settling 
chamber,  Bertoli  brand  extra  virgin  olive  oil  from  a  TSI  Six  Jet  Atomizer  Model  9306  feeds  a 
custom  built  seeder  rake  designed  and  constructed  by  Luker  (13)  and  Hale  (9)  and  modified  for 
this  thesis  (See  Figure  3.3).  According  to  the  Atomizer  Manual  (1),  the  average  particle  size  for 
olive  oil  is  0.6  pm. 


Settling  Chamber 


Figure  3.3  Seeder  System  Schematic 


The  air,  carrying  the  olive  oil  seed,  is  accelerated  through  the  nozzle  and  into  two  square  test 
sections  each  measuring  6.35  cm  by  6.35  cm  by  33.02  cm.  The  nozzle,  test  section,  and  diffuser  are 
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each  constructed  from  1.905cm  thick  aluminum  alloy.  The  test  section  is  modular  in  design  (Figure 
3.4),  allowing  different  components  to  be  mixed  and  matched.  Two  different  side  walls  were  used 
in  this  research.  To  allow  for  photographic  analysis  of  the  flow  field,  a  sidewall  was  constructed 
from  1.905  cm  thick  optical  grade  Plexiglass.  An  alternate  sidewall,  used  with  the  Laser  Doppler 
Velocimeter,  consists  of  an  aluminum  housing  protecting  two  3  inch  diameter  optical  pure  glass 
windows.  The  test  section  floors,  ceilings,  and  sidewalls  are  sealed  with  rubber  gaskets  to  prevent 
leakage.  Room  temperature  vulcanizer  was  also  used  for  any  remaining  test  section  leakage. 


Figure  3.4  Photograph  of  Test  Section  in  AFIT  Mach  3  Wind  Tunnel 


The  tunnel  air  exhausts  through  a  Kinney  Corporation  GP-8  gate  valve  into  the  vacuum 
system  with  a  tank  capacity  of  20  cubic  meters,  providing  an  approximate  run  time  of  35  seconds. 
The  tunnel  takes  approximately  8  seconds  to  reach  steady  state,  providing  usable  steady-state  flow 
for  roughly  25  seconds.  The  vacuum  is  provided  by  three  Model  212-11  MicoVac  vacuum  pumps 
manufactured  by  Stokes  Penwalt,  Figure  3.5.  Each  vacuum  pump  is  powered  by  a  7.6  metric 
horsepower  Reliance  electric  motor. 
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Figure  3.5  Vacuum  Pumps  for  AFIT  Mach  3  Wind  Tunnel 


3.3  Laser  Doppler  Velocimeter 

The  Laser  Doppler  Velocimeter  used  in  this  research  is  a  Dantec  Laser  Doppler  Anemometer 
with  the  Dantec  60X  series  fiber  flow  optics  system,  shown  in  Figure  3.6.  The  optical  system 
consists  of  a  laser,  routing  transmitter,  transmitting  optics,  receiving  optics  and  photomultiplier 
tubes.  The  system  uses  a  Ion  Laser  Technology  Model  5500A  laser  with  a  maximum  of  300  mW 
of  power  and  operates  in  the  frequency  range  from  457  nanometers  to  514.5  nanometers.  Luker 
(13)  showed  the  laser  has  degraded  in  efficiency  and  has  a  maximum  output  power  of  275mW.  This 
laser  is  powered  by  a  210V  and  20A  Laser  Technology  Model  5405A  power  supply  running  off  220 
VAC  power.  The  laser  beam  diameter  was  0.82  mm  and  the  laser  was  operated  in  the  TEMO0 
mode,  allowing  the  beam  to  be  centered  around  the  optical  axis  with  maximum  power  in  the  center 
of  the  beam. 

A  Fiber  Flow  60X41  transmitter  was  used  with  a  40  Mhz  Bragg  cell  that  split  the  laser  into 
514.5  nm  (green),  488.0  nm  (blue),  and  476.5  nm  (violet).  Each  wavelength  was  also  split  into 
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two  independent  beams  and  the  Bragg  cell  supplied  the  necessary  40  Mhz  frequency  shift  as  noted 
in  Luker  (13).  The  use  of  DANTEC  60X24  fiber  manipulators  allowed  the  fiberoptic  cables  to  be 
precisely  aligned  with  the  beam,  ensuring  maximum  power  and  efficiency.  These  cables  carried  the 
beam  to  the  transmitting  and  recieving  optics. 

The  collection  optical  probe  was  a  DANTEC  model  60X60  and  the  transmitting  probe  was 
a  DANTEC  model  60X61.  Each  probe  contains  transmitting  and  collecting  fibers,  an  optical 
transducer,  and  a  distribution  unit.  The  model  60X61  transmitting  probe  was  a  two-component, 
four-beam  probe  using  488  nm  and  514.5  nm  beams.  The  model  60X60  was  a  one-component, 
two-beam,  476.5  nm  probe.  The  transmitting  optics  were  mounted  on  the  traverse  mounting  bench 
using  DANTEC  Model  60X3611  probe  support  that  allowed  gross  angular  adjustments  and  the 
receiving  optics  were  mounted  on  a  DANTEC  Model  60X3631  probe  support  with  micrometer 
adjustments  for  precise,  three-dimensional  alignment. 

The  control  volume  generated  by  the  laser  beams  was  9.0  mm  across  the  tunnel  and  0.276 
mm  in  diameter.  Both  probes  used  a  600  mm  focal  length  lens  and  the  collection  optics  also 
used  a  DANTEC  Model  55X12  beam  expander,  allowing  additional  reflected  light  to  be  collected. 
The  collected  light  was  channeled  through  the  fiberoptic  cables  to  a  DANTEC  Model  55X35  color 
separator.  This  device  split  incoming  laser  energy  into  the  514.5  nm  and  the  488.0  nm  wavelengths 
and  routed  each  color  to  a  DANTEC  Model  55X35  Photomultiplier  tube.  These  tubes  converted 
the  incoming  light  into  electrical  signals  for  processing  by  a  Burst  Spectrum  Analyzer  (BSA). 

The  DANTEC  3D  Traverse  System  held  the  optics  and  was  used  to  position  the  lasers  for 
precise  control.  The  traverse  was  powered  by  three  Cleveland  Machine  Controls  stepper  motors 
regulated  by  Dynapar  encoders  Model  M20100003331  that  converted  1000  pulses  into  a  2mm  dis¬ 
placement  of  the  traverse  through  one  revolution  of  the  screw.  The  entire  system  accuracy  was 
rated  at  ±80  pm  over  600  mm  range  of  travel  and  could  support  a  maximum  load  of  445  N.  A  1.83m 


3-5 


DANTEC  laser  mounting  bench  replaced  the  traverse  mounting  bench  to  obtain  more  favorable 
optical  characteristics.  This  can  be  seen  in  figure  3.7. 


Figure  3.6  Laser  and  Photomultiplier  Tubes 


The  data  are  processed  by  a  Dantec  Burst  Spectrum  Analyzer  Model  57N20  Enhanced  (u 
component  of  velocity)  for  the  514. 5nm  (green)  laser  and  a  Dantec  Burst  Spectrum  Analyzer  Model 
57N35  (v  component  of  velocity)  for  the  488.0  nm  (violet)  laser.  The  57N20  is  the  master  BSA  for 
coincidence  filtering.  The  data  collection  is  performed  on  a  Gateway  2000  486DX/33Mhz  computer 
with  a  92032  IEEE  data  collection  card  feeding  the  Burst.ware  Version  3.00  data  acquisition  software 
system.  The  Burst  Spectrum  Analyzers  and  data  acquisition  computer  are  shown  in  figure  3.8. 

For  an  excellent  description  of  the  principles  of  Laser  Doppler  Velocimetry,  the  reader  is 
encouraged  to  see  Chapter  3  of  Luker  (13)  or  Chapter  3  of  Hale  (9). 
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3.4  Photographic  Analysis  Apparatus 


A  Xenon  Corporation  437B  Nanopulser  fed  a  DeLeone  N-787B  Broad  Spectrum  spark  gen¬ 
erator  providing  light  to  obtain  the  shadowgraphs  (See  Figure  3.9)  by  reflecting  the  light  from  a 
152.4  cm  focal  length/14.4  cm  diameter  mirror  and  using  type  57  ISO  3000/36°  Polaroid  film  and 
a  Polaroid  545i  camera  to  capture  the  image.  Using  an  additional  mirror  and  a  light-splitting  knife 
edge,  the  schlieren  images  were  obtained  using  the  same  light  source  and  an  additional  11.2  cm 
diameter  focusing  lens.  This  is  shown  in  Figure  3.10. 

Flow  Field 


Light  Source 


3.5  Pressure  and  Temperature  Measurement 

The  total  pressure  measurements  were  taken  with  an  Endevco  brand  0-0.69  MPa  pressure 
transducer  set  at  2.13  X  105±  690  Pa.  The  total  temperature  was  measured  using  an  Omega 
Engineering  Type  K  thermocouple. 

3.6  Model  Description 

The  model  used  to  generate  the  shock-boundary  layer  interaction  was  a  machined  aluminum 
wedge  rigidly  mounted  to  the  wind  tunnel  floor.  The  flow  deflection  angles  for  the  three  wedges 
were  5°,  7°,  and  10°.  Each  wedge  measured  7.21  cm  in  the  streamwise  direction  and  6.35  cm  across 
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Light  Source 


Flow  Field 


Figure  3.10  Schlieren  Setup 


the  wind  tunnel  test  section.  The  shock-generating  wedges  were  set  at  6.99  cm  downstream  of  the 
upstream  test  section  edge  so  that  the  shock  wave/boundary  layer  interaction  would  appear  in  the 
available  LDV  windows.  The  wedges  were  secured  to  the  tunnel  floor  with  four  screws  and  sealed 
with  a  rubber  gasket  as  seen  in  Figure  3.11.  This  designed  allowed  rapid  changeout  of  the  wedge 
without  removing  the  test  section  from  the  wind  tunnel. 


Figure  3.11  Photograph  of  Shock-Generating  Wedge  and  Tunnel  Floor 


Each  wedge  had  three  pressure  taps  equally  spaced  on  the  upstream  face.  The  tunnel  floor 
also  had  a  pressure  tap  1.27  cm  upstream  of  the  wedge  leading  edge.  The  tunnel  ceiling  had  nine 
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pressure  taps  spaced  1.0  cm  apart,  beginning  at  16.18  cm  from  the  upstream  edge  of  the  test 
section,  to  locate  the  shock  wave  reflection  from  the  ceiling.  The  pressure  ports  were  also  designed 
to  address  the  pressure  gradient  variations  between  shock  wave  angles  (See  Figure  3.12). 
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Pressure  Tap 


Incoming  Flow 


Incoming  How 


Wedge 


Incoming  How 


Figure  3.12  Pressure  Tap  Locations 


3.1  Wind  Tunnel  Test  Procedure 

When  the  vacuum  tanks  were  below  10  mmHg  the  tunnel  was  ready  for  operation.  The  first 
step  was  to  release  the  gate  valve  and  expose  the  test  section  to  the  vacuum  system.  Then,  the 
high  pressure  air  was  released  through  the  ball  valve  to  begin  high-speed  airflow.  The  tunnel  took 
approximately  8  seconds  to  reach  steady-state.  During  this  time,  two  jets  of  the  six-jet  seeder  system 
were  activated,  injecting  the  seed  into  the  test  section.  Once  steady-state  flow  was  achieved,  the 
Burstware  software  was  activated  to  collect  data  for  12  seconds,  sampling  as  much  data  as  possible. 
After  this  was  complete,  the  tunnel  was  shut  down  by  first  closing  the  ball  valve,  stopping  the  high 
pressure  air.  It  was  then  safe  to  close  the  gate  valve  and  begin  to  drawdown  the  vacuum  system 
again.  The  vacuum  system  could  recharge  in  approximately  6  minutes. 
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3.8  Traverse  and  Laser  Alignment 

To  simplify  LDV  operation  and  coordinate  transforms,  the  traverse  bench  was  set  perpendic¬ 
ular  to  the  test  section.  The  alignment  was  by  accomplished  by  aligning  the  lasers  parallel  to  the 
traverse  bench.  The  probe  supports  were  marked  at  1°  increments,  allowing  accurate  alignment 
with  the  bench.  The  traverse  assembly  was  then  rotated  until  the  lasers  reflected  back  on  them¬ 
selves,  ensuring  an  accurate  alignment.  Then,  the  lasers  were  rotated  3.5°  off  axis  of  the  bench, 
allowing  the  lasers  to  point  downstream  and  operate  in  the  forward-scatter  mode.  This  also  pro¬ 
tected  the  sensitive  photomultiplier  tubes  from  direct  laser  light  from  the  transmitting  lasers.  This 
alignment  does  cause  a  small  bias  in  the  data  but  it  was  considered  small  enough  to  be  neglected 
(See  Figure  3.13).  Notice  that  <f>  is  7°.  Luker  shows  that  a  7°  angle  provides  the  maximum  laser 
scattering  intensity(13). 


Test  Section 


Incoming  Flow 


Figure  3.13  Laser  and  Traverse  Alignment  Relative  to  Test  Section 

Since  the  shock  wave/boundary  layer  interaction  occurred  on  a  flat  plate,  the  problems  of 
LDV  angular  bias  were  minimized.  However,  to  place  the  origin  of  the  traverse  it  was  necessary 
to  align  the  lasers  with  the  test  section  floor.  Using  a  set  of  Argon-ion  laser  protective  goggles, 
the  laser  reflection  points  appeared  as  small  dots.  These  small  dots  were  located  at  the  point  of 
maximum  laser  intensity.  These  dots  were  placed  by  moving  the  traverse  such  that  the  two  514.5 
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nm  (green)  beams  were  bisected  by  the  tunnel  floor.  This  location  was  set  as  the  traverse  origin 
point. 

With  the  transmitting  optics  aligned  with  the  wind  tunnel,  the  receiving  optics  were  aligned 
with  the  control  volume.  Using  the  lasers  from  the  transmitting  optics,  the  two  control  volumes 
could  be  aligned  with  a  rough  degree  of  accuracy.  The  lasers  were  then  raised  to  the  freestream  of 
the  test  section  and  fine-tuned  by  injecting  a  small  amount  of  seed  into  the  test  section.  With  the 
tunnel  not  operating  and  using  the  on-line  display  feature  of  the  Burstware  software,  it  was  possible 
to  complete  the  alignment  process.  The  on-line  display  feature  constantly  updates  the  validity  and 
the  data  acquisition  rates  as  the  receiving  optics  were  adjusted.  This  allowed  the  receiving  optics 
to  be  adjusted  very  accurately  for  the  maximum  data  acquisition  and  validity  rates.  The  goal  for 
this  process  was  both  BSA’s  over  90%  validity  and  data  acquisition  rates  over  24  KHz.  It  was 
important  to  maximize  both  BSA’s  simultaneously  as  it  was  possible  to  achieve  high  rates  on  one 
BSA  while  the  other  BSA’s  rates  were  very  low. 


3.9  Burst  Spectrum  Analyzer  Settings 

The  fixed  software  settings  for  the  Burstware  Data  Acquisition  Software  are  presented  in 
Tables  3.1  to  3.3.  These  settings  were  kept  constant  throughout  the  experiment. 


Table  3.1  Quick  Menu  Constant  Settings 


Category 

BSA  1  Setting 

BSA  2  Setting 

Bandwidth  (m/s) 

325.12 

246.7 

Shifter  Mode 

Norm 

Norm 

Signal  Gain  (dB) 

42 

42 

High  Voltage  (V) 

1304 

1504 

Max.  Anode  Current  (mA) 

1.6 

1.6 

Pedestal  Attenuation  (dB) 

6 

6 

Duty  Cycle  (%) 

100 

100 

Dead  Time 

0 

0 
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Table  3.2  Soft  Menu  Constant  Settings 


Category 

BSA  1  Setting 

BSA  2  Setting 

Timer  Clock 

Master 

Slave 

Coincidence  Mode 

Master 

Slave 

Arrival  Time  Base 

Internal 

Internal 

Burst  Detection 

Both 

Both 

Oversize  Rejection 

1 

1 

Buffer  Mode 

FIFO 

FIFO 

Max.  Anode  Current  (mA) 

1.6 

1.6 

Quality  Factor  (%) 

50 

50 

Collection  Mode 

Burst 

Burst 

Table  3.3  BSA  Program  Menu  Constant  Settings 


Category 

BSA  1  Setting 

BSA  2  Setting 

Timeout  (Sec) 

12.0 

12.0 

Number  of  Bursts 

150000 

150000 

Velocity 

Yes 

Yes 

Transit  Time 

Yes 

Yes 

Arrival  Time 

Yes 

Yes 

In  addition  to  the  constant  settings  above,  the  record  length  was  kept  constant  at  16  samples. 
The  variable  settings  adjusted  during  the  experiment  were  the  center  frequency  on  BSA  1  (u- 
component  of  velocity),  and  the  size  of  the  FIFO  Buffer.  The  center  frequency  was  kept  as  close 
as  possible  to  the  average  velocity  at  any  given  point.  This  value  needed  to  be  adjusted  as  the 
point  of  interest  moved  from  the  wall  to  the  freestream.  The  size  of  the  FIFO  buffer  determined 
how  many  bursts  were  stored  on  the  BSA’s  before  transfer  to  memory.  If  the  buffer  size  was  set  at 
2000,  the  BSA’s  would  collect  2000  bursts  before  sending  them  to  memory.  The  problem  with  this 
buffer  was  the  last  group  of  bursts  would  not  be  sent  to  memory.  If  2500  bursts  were  collected,  the 
last  500  bursts  would  be  lost.  In  order  to  keep  as  many  bursts  as  possible,  the  FIFO  buffer  needed 
to  be  reduced.  However,  a  balance  was  sought  because  if  the  FIFO  buffer  was  too  small  the  BSA’s 
spent  too  much  time  transfering  data  and  not  collecting  it. 

Coincidence  filtering  was  performed  on  the  data  by  the  BSA’s  based  on  a  user-input  value. 
According  to  the  Burstware  Manual  (5),  best  results  are  achieved  when  the  coincidence  window  is 
set  at  twice  the  record  interval.  Realizing  that  at  650  m/s,  a  particle  will  transit  the  0.276mm  long 


3-14 


control  volume  in  roughly  0.4  ps,  following  this  advice  would  limit  the  record  length  to  8  samples. 
When  the  record  length  was  set  to  16  samples,  the  record  interval  was  0.333ps,  slightly  below  the 
minimum  transit  time.  The  effects  of  record  length  are  discussed  in  Luker  where  it  was  shown  that 
the  16  sample  record  length  was  well-behaved  while  the  8  sample  record  length  showed  enormous 
scatter  (13).  Therefore,  for  this  experiment,  the  coincidence  window  value  was  set  to  0.333ps  and 
the  record  length  was  set  to  16  samples. 
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IV.  Data  Reduction  Techniques 


4-1  Velocity  Calculation 

The  Dantec  Burst  Spectrum  Analyzer  used  a  Fast  Fourier  Transform  (FFT)  to  determine  the 
Doppler  frequency  of  the  scattered  laser  light.  This  frequency  was  converted  to  velocity  using  the 
frequency-velocity  calibration  factor,  Cjv.  During  this  project,  the  calibration  factors  for  the  u  and 
v  components  were 


Cfv,u  —  8-128^.^,  CfV)V 


7.709- 


MHz 


(4.1) 


The  total  velocity  can  be  shown  to  be  (13) 


Q=\fu- 


2  +  v2  +  (i>')2 


(4.2) 


and  the  Root  Mean  Square  (RMS)  value  of  the  velocity  components  becomes,  with  the  as¬ 
sumption  of  =  ~  1  and  =  <  1, 


(QO2  =  («')2  (4.3) 

4-2  Separation  of  Mean  and  Turbulent  Mach  Numbers 

As  previously  noted,  the  turbulent  characteristics  of  the  flow  greatly  complicate  the  calcula¬ 
tion  of  the  flow  variables.  The  Mach  number  is  no  different.  Generally,  the  Mach  number  is  given 
by 


a 


Q 


(4.4) 
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Expanding  the  Mach  number  into  mean  and  fluctuating  components  and  solving  for  M'  gives, 


after  manipulation  (13),  the  fluctuating  Mach  number  RMS  value 


TrWT 


(M')2  =  M^U- 
Q 


1  +  -(T-1  )M2 


(4.5) 


4-3  Mach  Number  Calculation 

The  local  Mach  number  was  calculated  using  the  measured  LDV  velocity  and  the  total  tem¬ 
perature  measured  in  the  plenum  along  with  the  assumptions  of  adiabatic  flow  and  a  calorically 
perfect  gas.  The  adiabatic  energy  equation  was  developed  for  turbulent  flow  by  expanding  the 
adiabatic  relationship  between  temperature,  stagnation  temperature,  and  Mach  number  (13).  The 
final  mean  Mach  number  expression  is 


M  = 


2  Q2 

2yRT0  -  Q2(7  -  1) 


(4.6) 


It  is  important  to  remember  that  this  Mach  number  expression  was  developed  assuming 
adiabatic  flow  (constant  total  temperature)  throughout  the  boundary  layer.  Kistler  (12)  verified 
this  assumption  by  measurement  on  a  flat  plate. 


4-4  Density  and  Temperature  Calculation 

Calculation  of  the  density  and  temperature  profiles  were  complicated  by  the  fact  that  the  only 
measurements  taken  in  the  test  section  were  the  wall  pressure  and  temperature.  The  density  and 
temperature  in  the  boundary  layer  were  estimated  based  on  certain  assumptions.  If  the  validity 
of  these  assumptions  breaks  down,  the  reported  values  of  the  temperature,  density,  and,  of  course, 
incompressible  Reynold’s  shear  stress  will  be  incorrect. 
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4-4-1  Flat  Plate.  On  the  flat  plate,  the  density  was  calculated  assuming  that  the  pres¬ 
sure  normal  to  the  surface  was  constant.  Using  this  assumption,  the  modified  Crocco-Busemann 
approximation  produces  (6) 


Pw 

P 


TV] 


=  l  +  B' 


(4.7) 


where  pw  and  Tw  are  the  density  and  temperature  at  the  wall  and  A'  and  B'  are  given  by 


{A')2  =  \  e 


1  +  ^ rM 2 
B  = - £ - -  -  1 


T*. 

Tc 


(4.8) 


(4.9) 


The  density  at  the  wall,  pw,  was  calculated  using  the  measured  wall  temperature  and  pressure 
and  the  perfect  gas  relation: 


_  Pw 
P w  pm 

ll'U  -*■  W 


(4.10) 


4-5  Boundary  Layer,  Displacement  and  Momentum  Thickness  Calculations 

The  boundary  layer  thicknesses,  based  on  velocity  (<5U),  was  calculated  for  the  flat  plate.  The 
boundary  layer  thicknesses  were  calculated  by  choosing  the  largest  velocity  from  a  given  traverse 
as  the  edge  value,  non-dimensionalizing  the  profile  by  the  edge  value,  and  searching  outward  from 
the  wall  until  the  value  of  ^  >  0.995  was  attained.  The  value  of  6  was  then  calculated  by 
linear  interpolation  between  the  data  points  just  above  and  below  the  boundary  layer  edge.  The 
displacement  and  momentum  thicknesses  were  calculated  using  the  equations  in  White  (19)  and 
implemented  by  Luker  (13). 
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4-6  Turbulence  Statistics 


Turbulence  statistics  are  needed  for  two  reasons.  First,  statistical  variables  appear  directly 
in  the  time-averaged  form  of  the  Navier-Stokes  equations.  In  addition,  the  turbulence  statistics 
help  to  understand  features  of  the  fluctuating  flowfield.  For  an  additional  discussion  of  turbulence 
statistics,  see  Luker  (13). 

The  fluctuation  from  the  mean,  x'  =  x—x,  is  the  variable  around  which  the  turbulent  statistics 
are  built.  Here,  the  variable  x  was  used  to  show  that  turbulent  statistics  were  generated  for  both 
the  u  and  v  components  of  velocity.  One  problem  appears,  though.  By  definition,  the  average  of 
x'  is  zero.  Therefore,  in  order  to  get  an  estimate  of  the  average  size  of  the  fluctuations,  the  Root 
Mean  Square  (RMS)  value  of  x'  was  used.  The  RMS  value  of  x'  is  defined  as 


Xrms  =  <7x  =  \J^T)  (4-n) 

Here  the  summation  is  performed  using  the  n  samples  collected.  Notice  that  the  denominator 
is  n  —  1  samples.  The  Burstware  software,  however,  calculates  the  RMS  value  over  n  samples.  See 
Appendix  B  for  a  more  complete  discussion  of  the  effect  of  this  error  on  the  results. 

The  Burstware  software  also  performs  calculations  to  determine  the  skewness  ( Skx )  and 
flatness  ( Flx )•  These  are  shown  below. 


Skx  = 


E(*')3 


(4.12) 


Fir  =  (4.13) 

The  skewness  is  important  when  using  LDV  because  it  helps  determine  if  the  results  are 
affected  by  seeding  problems.  On  the  other  hand,  the  flatness  is  used  to  calculate  the  intermittency 
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factor,  7U.  This  factor  is  the  ratio  of  the  flatness  of  the  Gaussian  curve  (3.0)  and  the  flatness  of 
the  measurements: 


*  =  K  <414> 

The  most  important  statistical  variable  is  the  velocity  correlation  u'v' .  This  appears  directly 
in  the  Reynold’s  Averaged  Navier-Stokes  equations  and,  along  with  density,  are  the  building  blocks 
of  the  incompressible  Reynolds  shear  stress.  The  correlation  was  calculated  using 


u'v' 


E(«v) 

n 


(4.15) 


Notice  that  while  the  averages  of  u'  and  v'  are  zero,  the  average  of  the  product  u'v'  may  not 
be  zero. 

The  final  statistical  value  calculated  was  the  velocity  correlation  coefficient,  Ruv,  which  is 
calculated  by 


Ruv  =  -  (4.16) 

RUv  was  used  to  provide  insight  in  the  nature  of  the  turbulence.  As  an  example,  if  the 
turbulence  was  isotropic  (no  preferred  direction  of  turbulent  structures),  Ruv  shold  have  had  a  value 
of  1.0.  Usually,  for  low  speed  flat  plate  flow  shows  that  the  magnitude  of  Ruv  is  approximately  0.4 
to  0.5.  (3) 

4-7  Wall  Shear  Stress 

The  wall  shear  stress  was  interatively  calculated  using  the  bisection  method  to  solve  the 
skin-friction  equation  developed  by  Van  Driest  (6)  and  implemented  by  Luker  (13). 
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where  k  =  0.41  and  lo  —  0.68.  It  should  be  noted  that  this  equation  was  developed  for 
turbulent  flow  over  a  flat  plate.  Rex  is  the  Reynolds  number  based  on  the  distance  from  the 
leading  edge  of  a  flat  plate  immersed  in  a  uniform  flowfield.  The  wind  tunnel  is  modeled  as  a  flat 
plate  with  the  leading  edge  at  the  nozzle  throat. 


For  a  discussion  of  the  Van  Driest  velocity  profile  as  implemented  in  this  study,  see  Chapter 
4  of  Luker  (13). 


4-8  Turbulent  Kinetic  Energy 

The  Turbulent  Kinetic  Energy  (TKE)  is  a  measure  of  the  energy  contained  in  the  turbulent 
motion  of  a  flowfield.  The  TKE  is  defined  as 


TKE  =  (u')2  +  (v')2  +  (w'f 
2 


(4.18) 


Even  in  a  two-dimensional  flowfield,  where  w  =  0,  the  lateral  fluctuation,  w',  is  not  necessarily 
zero.  With  the  current  experimental  setup  it  was  not  possible  to  measure  w' .  Therefore,  in  order 
to  calculate  TKE,  it  was  assumed  that  the  size  of  w'  is  the  same  order  as  v' ,  and  the  estimated 
TKE  relationship  becomes 


tke  =  ^1±^L ±oo! 
2 


(«') 


/\2 


A2 


+  (0 


(4.19) 
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V.  Results  and  Discussion 


5.1  Roadmap 

This  discussion  of  the  results  for  this  study  will  begin  with  an  explanation  of  the  coordinate 
systems  used  followed  by  a  summary  of  the  test  conditions.  Schlieren  and  shadowgraph  photographs 
for  each  flow  field  are  then  presented  in  the  section  on  flow  visualization.  The  flat  plate  Laser 
Doppler  Velocimetry  (LDV)  results  are  presented  and  compared  to  Luker  (13),  Hale  (9),  and 
Elena  and  LaCharme  (7).  The  shock  wave/boundary  layer  interaction  discussion  begins  with  a 
comparison  to  the  data  from  Rose  and  Johnson  (16)  for  a  similar  shock  wave/boundary  layer 
interaction  flowfield.  The  present  study  results  for  the  flowfield  induced  by  a  wall-mounted  5° 
wedge  are  then  presented,  to  illustrate  the  effects  of  a  shock  wave/boundary  layer  interaction. 
Finally,  the  effects  of  shock  strength  on  the  properties  of  the  flowfield  are  examined. 

5.2  Comments  on  Coordinate  Systems 

Throughout  this  report,  two  different  coordinate  axes  are  used.  In  order  to  avoid  confusion 
during  the  discussion,  each  will  be  developed. 

The  wind  tunnel  coordinate  scheme  is  a  Cartesian  system  with  the  tunnel  nozzle  at  x  =  0. 
The  y-coordinate  is  pointed  vertically  with  the  positive  direction  down.  The  z-coordinate  completes 
the  axis  system  using  the  right-hand  rule. 

The  shock  wave/boundary  layer  interaction  coordinate  system,  denoted  by  the  capital  letters 
(X,  Y ,  Z ),  are  used  to  describe  variables  in  relation  to  the  shock  wave/boundary  layer  interaction. 
The  X-coordinate  is  pointed  down  the  tunnel  with  positive  in  the  downstream  direction.  The  Y- 
coordinate  is  pointed  positive  down  and  the  Z-coordinate  completes  the  system  using  the  right-hand 
rule  (see  Figure  5.1). 
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Figure  5.1  Wind  Tunnel  and  Shock  Wave/Boundary  Layer  Interaction  Coordinate  Systems 
5.3  Test  Conditions 

The  pressure  transducer  located  in  the  settling  chamber  showed  that  the  wind  tunnel  was 
operating  with  a  stagnation  pressure  of  2.18X105  ±  3.4X103  Pa.  The  stagnation  temperature 
measured  in  the  settling  chamber  was  297K  ±  2. IK.  This  variation  was  averaged  between  traverses. 
The  stagnation  temperature  climbed  throughout  a  single  traverse  but  the  variation  was  usually  less 
than  IK.  The  temperature  climb  was  attributed  to  the  operation  of  the  compressors.  The  mean 
temperature  never  rose  above  300K.  Because  of  this  temperature  rise,  the  stagnation  pressure  and 
temperature  were  recorded  during  each  data  point  and  used  during  data  reduction  accordingly.  The 
Mach  number  of  the  flow  was  approximately  2.9  with  Rex  —  8.423X106  and  Reso  =  1.332X105. 

5-4  Flow  Visualization 

Flow  visualization  photographs,  both  Schlieren  and  shadowgraphs,  were  taken  of  each  flowfield 
using  the  setup  described  in  Chapter  4.  The  flow  for  each  photo  is  left  to  right. 

5-4-3  Flat  Plate.  The  flat  plate  Schlieren  photograph  is  shown  in  Figure  5.2.  The  waves 
crossing  the  tunnel  in  an  ”  X”  pattern  are  shocks  caused  by  slight  misalignments  of  the  tunnel  test 
sections. 
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Figure  5.2  Schlieren  of  Flat  Plate  Using  Horizontal  Knife  Edge 
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The  flat  plate  shadowgraph  photo,  shown  in  Figure  5.3,  clearly  shows  that  the  seam  shocks 
are  very  weak  and  have  little  effect  on  the  flow. 


Figure  5.3  Shadowgraph  of  Flat  Plate  Flowfield 


5.5  Effective  Flow  Turning  Angle 

During  initial  model  checkout  runs  of  the  wind  tunnel,  it  was  observed  that  the  shock  wave 
interaction  was  occurring  significantly  upstream  of  the  location  expected  by  inviscid  shock  wave 
theory.  As  this  moved  the  interaction  outside  of  the  available  LDV  viewing  area,  this  was  of  great 
concern.  The  shock  wave  generators  were  designed  assuming  inviscid,  M  =  2.90,  uniform  flow. 
However,  from  initial  shadowgraph  photographs,  the  shock  wave  angle  for  the  5°,  7°,  and  10° 
were  26.5°,  28.5°,  and  31.5°  with  an  error  of  ±0.5°.  Initally,  the  Mach  Number  of  the  tunnel  was 
considered  the  suspect  value.  However,  after  examination  of  the  tunnel  and  comparing  the  present 
configuration  with  other  work  performed  in  this  facility,  the  Mach  Number  was  ruled  out.  After 
careful  consideration,  the  suspect  value  is  the  assumption  of  inviscid  flow.  It  appears  the  boundary 
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layer,  on  the  same  order  as  the  shock  generating  wedge,  causes  the  flow  to  turn  more  sharply  for 
a  given  wedge  angle.  Using  the  6-f3- Mach  relationship  (2)  and  holding  the  Mach  Number  at  2.90 
and  using  the  above  values  for  the  shock  wave  angle,  the  effective  turning  angles,  denoted  by  9eff, 
became  8.17°,  10.40°,  and  13.52°,  respectively.  It  appears  that  the  boundary  layer  constricts  the 
flow  into  a  smaller  area,  producing  a  greater  turning  angle.  Even  though  the  effective  turning  angle 
is  greater  than  the  wedge  angle,  for  clarity,  the  wedge  angle  will  be  continue  to  be  used  to  identify 
each  flowfield. 

5.5.1  Shock  Wave/Boundary  Layer  Interaction.  The  shock  wave/boundary  layer  interac¬ 
tion  for  the  shock  wave  induced  by  the  5°  wedge  is  shown  Figure  5.4.  The  5°  wedge  flow,  because 
of  the  low  shock  wave  angle,  needed  two  window  views  in  order  to  visualize  the  entire  interaction 
region.  Figure  5.4  shows  the  interaction  in  the  upper  right  corner  of  the  photo.  The  accompanying 
shadowgraph  is  shown  in  Figure  5.5.  The  shock  wave  in  the  lower  half  of  the  photos  is  the  shock 
wave  off  the  trailing  edge  of  the  wedge,  after  which  the  flow  is  again  aligned  with  the  tunnel  walls. 

The  downstream  section  of  the  5°  wedge  flow  is  shown  in  Figures  5.6  and  5.7.  The  interaction 
is  in  the  extreme  upper  left  corner  of  the  photos. 

The  shock  wave/boundary  layer  interaction  for  the  7°  wedge  flowfield  is  shown  in  Figures  5.8 
and  5.9.  The  larger  shock  wave  angle  of  the  7°  wedge  flowfield  have  allowed  the  entire  interaction 
region  to  be  viewed  in  one  photo.  Note  the  additional  waves  emanating  from  the  boundary  layer 
after  the  interaction  in  Figure  5.9.  Also,  note  the  clear  intermittency  at  the  top  of  the  boundary 
layer  at  the  lower  edge  of  Figure  5.8.  At  the  vertex  of  the  shock  wave  reflection,  note  the  small 
lambda  shock  developing.  From  this  photo,  it  can  be  seen  that  the  7°  wedge  flow  may  have  a  small 
separated  region  at  the  interaction.  However,  for  the  5°  wedge  flow  there  is  no  evidence  of  a  lambda 
shock  wave,  allowing  the  conclusion  to  be  drawn  that  the  5°  wedge  flow  is  not  separated  at  the 
interaction. 
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Figure  5.4  Schlieren  of  5°  Wedge  Flow  Interaction  Using  Horizontal  Knife  Edge  (Forward 
Window) 
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Figure  5.5  Shadowgraph  of  5°  Wedge  Flow  Interaction  Flowfield  (Forward  Window) 
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Figure  5.6  Schlieren  of  5°  Wedge  Flow  Interaction  Using  Horizontal  Knife  Edge  (Rear  Window) 
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Figure  5.7  Shadowgraph  of  5°  Wedge  Flow  Interaction  Flowfield  (Rear  Window) 
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Figure  5.9  Shadowgraph  of  1°  Wedge  Flow  Interaction  Flowfield 
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Figures  5.10  and  5.11  show  the  shock  wave/boundary  layer  interaction  for  the  10°  wedge 
flowfield.  Both  of  these  figures  show  a  small  lambda  shock  forming  at  the  interaction.  Again,  note 
the  additional  waves  emanating  from  the  boundary  layer  after  the  shock  wave  reflection.  In  these 
figures,  the  lambda  shock  pattern  at  the  vertex  of  the  shock  waves  is  sharper  and  larger.  The  10° 
wedge  flow  has  a  larger  separated  region  than  the  7°  wedge  flow. 


Figure  5.10  Schlieren  of  10°  Wedge  Flow  Interaction  Using  Horizontal  Knife  Edge 


5.6  Pressure  and  Temperature  Measurements 

The  pressure  distribution  over  the  flat  plate  with  the  impinging  shock  wave  is  shown  in 
Figure  5.12.  Notice  the  sharp  rise  in  pressure  where  the  shock  wave  impinges  on  the  flat  plate. 
The  flat  plate  pressure  distribution  is  also  presented  for  comparison.  Downstream  of  the  boundary 
layer/shock  wave  interaction,  the  pressure  again  drops  towards  the  freestream  value  after  the  shock 
wave  reflects  off  the  boundary  layer.  As  expected,  the  pressure  rise  increases  with  increasing  wedge 
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angle  and  is  greatest  for  the  10°  wedge  flowfield.  The  pressure  rise  from  the  shock  waves  generated 
by  the  1°  and  the  5°  wedges  are  less  severe  as  shown  in  Figure  5.12.  The  pressure  rise  normalized 
by  the  first  pressure  measurement  is  presented  in  Table  5.1. 


Figure  5.12  Pressure  Distribution  over  Flat  Plate  with  Impinging  Shock  Wave 


10° 

7° 

5° 

2.34 

2.02 

1.91 

Table  5.1  Normalized  Shock  Wave  Interaction  Pressure  Increase 

The  pressure  distribution  on  the  shock-generating  wedge  is  shown  in  Figure  5.13.  As  expected, 
the  10°  wedge  had  the  highest  pressure  followed,  respectively,  by  the  7°  and  the  5°  wedges. 

5.7  Flat  Plate  LDV  Results 

In  order  to  ensure  that  the  LDV  system  and  wind  tunnel  were  working  correctly,  it  was 
important  to  compare  to  earlier  studies.  The  works  of  Luker  (13)  and  Hale  (9)  were  chosen  for  a 
number  of  reasons.  The  author  assisted  on  the  data  collection  for  Luker  and  Hale  while  learning  the 
tunnel  and  LDV  systems,  and  the  data  were  readily  available.  Both  Luker  and  Hale  did  exhaustive 
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Figure  5.13  Pressure  Distribution  over  Shock  Generating  Wedges 


studies  of  the  flat  plate  characteristics  of  the  AFIT  Mach  3  Wind  Tunnel  and  this  data  are  used 
to  show  the  accuracy  of  the  present  system.  The  data  are  also  compared  to  Elena  and  LaCharme 
1988  flat  plate  experiments  (7).  In  addition,  the  flat  plate  flowfield  serves  as  a  baseline  against 
which  the  effects  of  the  shock  wave/boundary  layer  interaction  may  be  ascertained. 

5.7.1  Flat  Plate  Mean  Flow  Characteristics.  Mean  flow  properties  are  compared  to 
Luker’s  (13)  and  Hale’s  (9)  data  in  the  form  of  normalized  velocity  profiles  jj->  and  Mach  number 
profiles. 


5. 7. 1.1  Flat  Plate  Velocity  Profiles.  To  determine  the  flat  plate  boundary  layer 
thickness,  the  streamwise  mean  velocity  was  scaled  by  the  edge  value  (603.16  m/s  at  x  —  48  cm) 
and  plotted  versus  normalized  distance  from  the  wall  (Figure  5.14).  It  is  important  to  scale  the 
velocity  data  by  the  edge  value  to  eliminate  any  variations  in  test  conditions  between  test  traverses 
across  the  boundary  layer.  The  boundary  layer  thickness  was  taken  as  the  distance  from  the  wall 
where  jF  =  0.995.  The  flat  plate  boundary  layer  thickness  for  this  experiment  was  found  to  be 
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roughly  7.6  mm.  The  agreement  with  Luker’s  (13)  and  Hale’s  (9)  measurements,  seen  in  Figure 
5.14,  is  excellent. 


Figure  5.14  Normalized  Flat  Plate  Velocity  Profile 

5 .7. 1.2  Flat  Plate  Mach  Number.  The  Mach  number  was  calculated  as  described 
in  Chapter  4  and  is  shown  in  Figure  5.15.  The  data  presented  here  is  for  x  =  48  cm,  equating  to 
approximately  63.25  boundary  layer  thicknesses  (6jp)  downstream  of  the  wind  tunnel  nozzle  throat. 
The  comparison  data  from  Luker  and  Hale  was  taken  at  x  =  71  cm  (94.67  6fp  downstream  of  the 
nozzle  throat).  Luker  commented  on  the  Mach  number  decay  with  downstream  distance.  This  is 
seen  in  Figure  5.15.  The  present  study  is  23  cm  upstream  of  Luker’s  (13)  and  Hale’s  (9)  and  the 
Mach  number  is  0.1  higher. 

The  fluctuating  profile  comparison  is  presented  in  Figure  5.16.  In  general,  there  is  good  agree¬ 
ment.  The  present  data  tends  toward  a  higher  fluctuating  Mach  number  but  it  never  approaches  M 
=  1.0,  eliminating  any  tendency  to  produce  shocklets  (18).  In  addition,  Morkovin’s  hypothesis  can 
be  stated  as  ’’the  dynamics  of  a  compressible  boundary  layer  will  follow  the  incompressible  pattern 
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Figure  5.15  Flat  Plate  Mach  Number  Profile 

closely,  as  long  as  the  fluctuating  Mach  number,  M' ,  remains  small.”  (18)  Since  the  M'  for  the  flat 
plate  in  the  present  study  never  exceeds  0.35,  Morkovin’s  hypothesis  should  hold. 

5.7.2  Flat  Plate  Velocity  Fluctuations.  The  streamwise  velocity  fluctuations,  scaled  versus 
both  local  and  edge  velocities,  are  presented  in  Figures  5.17  and  5.18.  Note,  in  Figure  5.17,  the 
excellent  agreement  with  both  Luker  and  Hale  (13)  (9)and  Elena  and  LaCharme  (7).  Although  the 
comparison  to  Robinson,  et  al.,  (15)  was  not  as  good,  it  was  still  reasonable. 

The  cross-stream  velocity  fluctuation  is  presented  in  Figures  5.19  and  5.20,  again  scaled  by 
local  velocity  and  edge  velocity,  respectively.  When  scaled  against  the  local  velocity,  the  agreement 
appears  excellent.  However,  when  compared  to  the  data  of  Robinson,  et  al.,  (15),  the  agreement 
noticeably  worsens.  The  present  data,  along  with  Luker’s  (13)  and  Hale’s  (9)  cross-stream  velocity 
fluctuation,  appear  to  fall  much  lower  closer  to  the  wall.  However,  in  general,  the  agreement  with 
the  literature  is  considered  acceptable. 
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Figure  5.17  Flat  Plate  Streamwise  Velocity  Fluctuations  Scaled  by  Local  Velocity 
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Figure  5.18  Flat  Plate  Streamwise  Velocity  Fluctuations  Scaled  by  Edge  Velocity 
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Figure  5.19  Flat  Plate  Cross-stream  Velocity  Fluctuation  Scaled  by  Local  Velocity 
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Figure  5.20  Flat  Plate  Cross-stream  Velocity  Fluctuation  Scaled  by  Edge  Velocity 

5.7.3  ’’Incompressible”  Reynolds  Shear  Stress.  The  incompressible  Reynolds  shear  stress 
is  defined  as  —'pu'v' .  To  facilitate  comparison  with  other  studies,  it  is  important  to  represent  the 
shear  stress  in  a  non-dimensional  fashion.  By  noting  that  -^v'  is  equivalent  to  the  incompressible 
Reynold’s  Shear,  —'pu'v' ,  divided  by  pu2,  a  convenient  and  accurate  scaling  is  obtained.  This  scaling 
is  important  for  three  reasons.  First,  this  method  of  scaling  eliminates  any  linear  biasing  which 
may  influence  the  measurements  (see  Appendix  B  of  Luker  (13)).  Also,  this  scaling  allows  easier 
comparison  to  CFD  results,  because  turbulence  models  cannot  resolve  u'v'  but  can  resolve,  rjy  m 
—pu'v1.  More  importantly,  however,  this  scaling  allowed  the  presentation  of  the  incompressible 
turbulent  shear  stress  profile  using  only  directly  measured  information.  This  eliminates  errors 
occurring  from  inaccurate  calculations  of  other  scaling  properties.  Figure  5.21  shows  excellent 
agreement  between  Luker’s  (13)  and  Hale’s  (9)  data  and  the  present  study  for  the  ’’incompressible” 
shear  stress. 

The  flat  plate  results  presented  above  are  considered  acceptable  and  allow  consideration  of 
the  primary  flowfield  of  interest,  the  shock  wave/boundary  layer  interaction. 
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Figure  5.21  Incompressible  Reynolds  Shear  Stress  Scaled  by  Local  Velocity 


5.8  Shock  Wave/Boundary  Layer  Interaction  LDV  Results 

5.8.1  Scaling  Factors.  During  the  present  study  it  was  noted  that  the  boundary  layer 
thickness  at  each  location  would  be  extremely  difficult  to  obtain  for  a  number  of  reasons.  Most 
importantly,  the  shock  wave  itself  interferes  with  the  boundary  layer  at  the  upstream  measurement 
stations.  In  addition,  the  velocity  profiles  at  the  downstream  stations  did  not  have  a  classic  bound¬ 
ary  layer  profile.  Therefore,  it  was  decided  to  scale  the  data  by  a  universal  set  of  scale  factors.  In 
addition,  to  facilitate  comparison,  a  similar  set  of  scale  factors  was  developed  for  the  data  of  Rose 
and  Johnson  (16).  These  values  are  shown  in  Table  5.2. 


Present  80  (mm) 

Present  U0  (m/s) 

Me 

n  M. 

re  m3 

7.59 

608.0 

2.87 

0.224 

Rose  and  Johnson  80  (mm) 

Rose  and  Johnson  U0  (m/s) 

22.9 

613.0 

Table  5.2  Scaling  Factors  for  Present  Study  and  Rose  and  Johnson 


5.8.2  Comparison  to  Rose  and  Johnson  Shock  Wave/Boundary  Layer  Interaction.  The 
work  of  Rose  and  Johnson  (16)  was  chosen  for  comparison  because  of  the  experimental  similarity. 
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Their  experiment  turned  the  flow  7°,  creating  a  shock  wave  that  impinged  on  a  flat  plate,  and 
studied  the  turbulence  properties  both  upstream  and  downstream  of  the  interaction  using  both 
LDV  and  hot-wire  anemometry. 

As  seen  in  Figure  5.22,  there  is  good  agreement  with  the  velocity  profile  data  of  Rose  and 
Johnson  (16)  at  the  upstream  station.  See  Section  5.5  for  a  discussion  on  the  effective  turning 
angle,  9ejf  ■  It  is  important  to  note  that  Rose  and  Johnson’s  upstream  measurement  station  was  at 
13.7  cm  or  X  =  1.6 660  upstream  of  the  shock  wave/boundary  layer  interaction  while  the  present 
study’s  upstream  location  was  1.0  cm  or  X  —  1.32 80  upstream  of  the  vertex  of  the  shock  wave 
reflection. 


Figure  5.22  Normalized  Velocity  Profile  Comparison  at  Upstream  Stations 


1st  Station 

2nd  Station 

1  st  Rose  and  Johnson  Station 

2nd  Rose  and  Johnson  Station 

1.32  60 

2.646„ 

1.6660 

2.75<50 

Table  5.3  Downstream  Measurement  Stations  in  Boundary  Layer  Thicknesses  from  Interaction 


The  velocity  profile  agreement  between  the  present  study  and  Rose  and  Johnson  experiment 
is  good.  The  measurements  from  Rose  and  Johnson  at  21.3  cm  or  X  =  1.6660  downstream  of  the 
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The  agreement  for  the  streamwise  velocity  fluctuations  at  the  upstream  location  between  the 
present  study  and  Rose  and  Johnson  is  fair.  Rose  and  Johnson’s  streamwise  velocity  fluctuation 
data  is  consistently  lower  than  the  present  study  but  does  follow  the  same  trend.  In  Figure  5.25, 
note  the  three  data  points  of  Rose  and  Johnson  near  the  value  of  j-  =  1.0.  At  similar  locations,  all 
the  data  traces  are  abruptly  shifted  towards  zero.  This  is  an  indication  of  shock  passage  through 
the  data  traces. 


Figure  5.25  Streamwise  Velocity  Fluctuation  Comparison  at  the  Upstream  Station 


On  the  other  hand,  the  downstream  streamwise  turbulence  intensity  profiles,  in  Figures  5.26 
and  5.27,  show  excellent  agreement  within  an  acceptable  scatter.  Both  studies  have  similar  mea¬ 
surement  distances  downstream  of  the  shock  wave/boundary  layer  interaction  and  this  fact  bears 
out  in  the  data. 

The  cross-stream  velocity  fluctuations  at  the  upstream  measurement  stations,  shown  in  Figure 
5.28,  also  shows  excellent  agreement.  Rose  and  Johnson’s  data  follows  the  pattern  for  the  7°  wedge 
closely  and  all  data  points  from  the  present  study  are  clustered  around  Rose  and  Johnson’s  data. 
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Figure  5.26  Streamwise  Velocity  Fluctuation  Comparison  at  ls<  Downstream  Stations 
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Figure  5.27  Streamwise  Velocity  Fluctuation  Comparison  at  the  2nd  Downstream  Stations 


Figure  5.28  Cross-stream  Velocity  Fluctuations  Comparison  at  the  Upstream  Station 


On  the  other  hand,  the  downstream  agreement,  shown  in  Figures  5.29  and  5.30,  is  only  fair. 
The  present  study’s  magnitude  is  greater  but  roughly  equal  in  scatter. 

For  the  Reynolds  shear  stress  profiles,  Rose  and  Johnson’s  data  are  scaled  by  the  wall  shear 
stress  value  reported  in  reference  (16).  These  values  are  157  ^  for  the  X  —  —  1 .66<50  station  and 
221  at  the  X  =  2.75<50  station.  The  Reynolds  shear  stress  upstream  station  comparison  is  shown 
in  Figure  5.31.  The  agreement  at  the  upstream  stations  is  excellent  between  Rose  and  Johnson’s 
data  and  the  present  study.  The  calculation  of  tw  used  for  this  comparison  is  discussed  later  in 
this  chapter  (Section  5.9.2). 

However,  the  Reynolds  shear  stress  agreement  at  the  downstream  stations  is  poor.  The 
present  data  for  Reynold’s  shear  stress  are  scaled  by  the  wall  shear  stress,  tw.  The  calculation 
of  tw  downstream  of  the  interaction  in  this  study  is  considered  suspect  and  is  a  poor  choice  as  a 
scaling  factor  but  Rose  and  Johnson  do  not  report  any  other  adequate  scaling  factor  in  (16),  choosing 
instead  to  leave  the  Reynolds  shear  stress  in  dimensional  form.  It  appears  that  the  assumption  of 
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Figure  5.29  Cross-stream  Velocity  Fluctuation  Comparison  at  Ist  Downstream  Station 
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Figure  5.30  Cross-stream  Velocity  Fluctuation  Comparison  at  2nd  Downstream  Station 
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_ Vt. _ 

Figure  5.31  Reynolds  Shear  Stress  Scaled  by  Wall  Shear  Stress  for  the  Upstream  Station 

Couette  flow  may  fail  at  this  measurement  station.  The  calculation  of  tw  is  explained  in  Section 
5.9.2. 

The  agreement  with  Rose  and  Johnson’s  (16)  data  is  good  and  the  studies  do  not  blatantly 
disagree  or  contradict  each  other,  with  the  possible  exception  of  the  downstream  Reynolds  shear 
stress  scaled  by  the  wall  shear  stress.  This  has  provided  a  sound  foundation  for  the  discussion  of 
the  effects  of  shock  strength  and  measurement  location  on  the  turbulence  properties  of  a  boundary 
layer  subject  to  a  shock  wave  interaction. 

5.8.3  Shock  Wave  Boundary  Layer  Results.  This  section  will  cover  the  shock  wave/boundary 
layer  results  for  the  5°  wedge  flow  to  illustrate  the  effects  of  a  shock  wave  impingement  on  a  turbu¬ 
lent  supersonic  boundary  layer.  Initially,  the  mean  flow  characteristics  are  discussed,  followed  by 
an  explanation  of  the  turbulent  properties. 
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Figure  5.32  Reynolds  Shear  Stress  Scaled  by  Wall  Shear  Stress  for  2nd  Downstream  Station 

The  velocity  profiles  were  obtained  by  dividing  the  streamwise  mean  velocity  by  the  velocity 
reference  value,  U0.  These  results  are  shown  plotted  versus  j-  to  facilitate  comparison  and  greatly 
simplify  the  data  reduction  process. 

The  effects  of  measurement  location  on  the  velocity  profiles  for  the  5°  wedge  flowfield  is  shown 
in  Figure  5.33.  The  upstream  station,  at  X  =  — 1.3260,  follows  the  flat  plate  data  closely.  The 
profiles  at  X  =  0<5o,  X  =  1.32<50,  and  X  =  2.64 50  match  each  other  closely,  showing  the  effect  of 
the  shock  wave  on  the  velocity  downstream.  It  appears  that  the  velocity  profiles  remain  relatively 
unchanged  downstream  of  the  interaction,  considering  each  profile  is  scaled  by  a  reference  quantity 
and  not  an  edge  value.  Notice  the  profile  at  X  =  —1.32  80  rapidly  slows  above  ¥-  ps  1.0,  where  the 
shock  wave  passed  through  the  traverse  profile. 

The  Mach  number  profiles  for  the  5°  wedge  flow  are  shown  in  Figure  5.34.  The  downstream 
Mach  numbers  are  roughly  0.2  lower  than  the  flat  plate  or  the  upstream  values.  The  stations  at 
0 60,  1.32<50,  and  2.64<50  are  closely  matched  and  show  similar  trends  to  the  velocity  profiles. 
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Figure  5.34  5°  Wedge  (6efj  =  8.2°)  Mach  Number  Profiles  at  each  Measurement  Station 


As  seen  in  Figure  5.35,  the  streamwise  velocity  fluctuations  are  increased  by  the  shock  wave 
interaction.  Near  the  wall,  the  increase  is  approximately  29%.  Further  from  the  wall  the  increase 
is  smaller  but  still  distinct.  There  is  also  some  upstream  influence  near  the  wall  where  the  X  = 
—  1.32^  station  streamwise  velocity  fluctuations  are  larger  than  the  flat  plate. 
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Figure  5.35  5°  Wedge  (9eff  =  8.2°)  Streamwise  Velocity  Fluctuations 

Like  the  streamwise  fluctuations,  the  cross-stream  velocity  fluctuations,  as  seen  in  Figure 

5.36,  are  energized  by  the  shock  wave  interaction  and  remain  elevated  through  the  X  =  2.64<50 
measurement  station.  However,  no  distinct  trend  based  on  downstream  location  can  be  discerned 
from  these  data.  Note  that  the  upstream  profile  is  intersected  by  the  shock  at  roughly  X  =  0.9 <5,, 
and  values  above  this  should  not  be  considered  when  comparing  the  upstream  to  downstream 
stations. 

Along  the  same  lines,  the  shock  wave/boundary  layer  interaction  has  a  distinct  effect  on  the 
Reynolds  shear  stress.  The  downstream  measurement  locations  have  a  much  higher  shear  stress 
near  the  wall  than  the  upstream  station.  The  increase  for  the  5°  wedge  flow,  as  shown  in  Figure 

5.37,  equates  to  nearly  a  56%  rise  in  shear  stress  downstream  of  the  interaction  over  the  upstream 
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Figure  5.36  5°  Wedge  (Oeff  =  8.2°)  Cross-stream  Velocity  Fluctuations 

values.  The  shear  stress  gradients  perpendicular  to  the  wall  are  much  higher  downstream  of  the 
interaction,  allowing  the  profiles  to  return  to  upstream  values  by  a  height  of  Y  =  O.5<50 . 

Concluding  the  examination  of  the  effect  of  the  shock  wave/boundary  layer  interaction  on  the 
5°  wedge  flow,  the  TKE  profiles  are  presented  in  Figure  5.38.  The  interaction  has  energized  the 
energy  in  the  flow  downstream,  increasing  the  TKE  at  all  downstream  locations.  The  TKE  near 
the  wall  increases  with  downstream  distance,  reaching  a  maximum  value  at  X  =  2.64<50,  the  final 
measurement  station. 

5.9  Effect  of  Shock  Strength 

To  examine  the  effect  of  shock  strength  on  the  turbulent  properties  of  a  shock  wave/boundary 
layer  interaction,  three  wedge  angles  were  used:  5°,  7°,  and  10°.  Streamwise  profiles  were  obtained 
at  a  constant  height  above  the  wall  of  Y  =  0.666o  will  be  used.  This  profile  orientation  facilitated 
the  analysis  of  shock  strength,  along  with  the  effect  of  downstream  distrance,  on  the  evolving 
boundary  layer. 
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The  velocity  defect  through  the  shock  wave/boundary  layer  interaction  is  dependent  on  shock 
strength,  as  seen  in  Figure  5.39.  The  velocity  drop  occurs  slightly  upstream  of  the  interaction  and 
continues  downstream  roughly  O.4<50.  Note  that  the  10°  wedge  profile  never  reachs  a  steady  state 
velocity,  continuing  to  accelerate  through  the  last  measurement  at  X  =  A.Q80. 


Figure  5.39  Streamwise  Velocity  Profile  at  Height  from  Wall  Y  =  0.66<5o 


Figure  5.40  shows  the  Mach  number  profile  versus  distance  from  the  shock  wave/boundary 
layer  interaction  at  a  height  from  the  wall  of  Y  =  0.666c,.  This  figure  clearly  shows  the  effect  of 
increasing  shock  strength  on  the  Mach  number  as  the  flow  moves  through  the  interaction.  The 
Mach  number  gradient  is  greater  for  increasing  shock  strength.  Figure  5.40  also  shows  that,  in  the 
10°  wedge  case,  the  Mach  number  has  recovered  from  the  interaction  and  is  approaching  a  new 
equilibrium  value  by  roughly  2.06o  downstream. 

The  streamwise  velocity  fluctuations  for  the  different  shock  strengths  at  a  height  from  wall 
of  0.666o  is  seen  in  Figure  5.41.  These  data  support  the  conclusion  that  increasing  shock  strength 
causes  an  increase  in  the  streamwise  velocity  fluctuations.  The  nature  of  the  10°  wedge  curve  also 
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Figure  5.40  Streamwise  Mach  Number  Profiles  at  Height  from  Wall  Y  =  0 ,6680 


supports  the  conclusion  that  this  flow  is  fundamentally  different  than  either  the  5°  or  the  7°  wedge 
flows.  The  evidence  supports  that  the  10°  flow  has  separated  at  the  interaction. 

Similar  to  the  streamwise  velocity  fluctuations  (Figure  5.42),  the  cross-stream  velocity  fluctu¬ 
ations  increase  as  a  function  of  shock  strength  near  the  interaction.  Farther  downstream,  however, 
this  conclusion  no  longer  seems  valid.  However,  the  downstream  values  for  cross-stream  veloc¬ 
ity  fluctuations  remain  elevated  over  the  upstream  values  through  the  last  measurement  location, 
X  =  4.0<5o. 

The  incompressible  Reynolds  shear  stress  is  defined  as  —pu'v1.  However,  to  facilitate  com¬ 
parison  with  other  studies,  it  is  important  to  represent  the  shear  stress  in  a  non-dimensional 
fashion.  The  correlation  is  presented  below  for  each  shock  strength  thoughout  the  shock 

wave/boundary  layer  interaction.  It  is  hoped  that  this  representation  of  the  incompressible  Reynolds 
shear  stress  will  allow  easy  comparison  to  CFD  results.  Other  reasons  for  this  scaling  are  outlined 
in  flat  plate  Reynolds  shear  stress  discussion  in  Section  5.7.3. 
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Figure  5.41  Streamwise  Velocity  Fluctuations  at  Height  from  Wall  Y  =  0.666, 


Figure  5.42  Cross-stream  Velocity  Fluctuations  at  Height  from  Wall  Y  =  0.665, 


Examining  the  incompressible  Reynolds  shear  stress  streamwise  profile  at  Y  =  0 ,66S0  (see 
Figure  5.43),  the  profiles  for  the  5°  and  the  7°  wedge  drop  just  before  the  interaction  and  return  to 
near  upstream  values  downstream  of  the  interaction.  The  10°  wedge  flow,  however,  rises  sharply 
just  before  the  interaction  and  drops  below  the  upstream  values  downstream  of  the  interaction. 
Thus,  it  appears  the  character  of  the  flow  may  have  changed  with  the  higher  shock  strength. 
The  10°  wedge  flow  may  have  separation  near  the  interaction,  as  suggested  by  the  Schlieren  and 
shadowgraph  photos  in  Figures  5.10  and  5.11. 
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Figure  5.43  Incompressible  Reynolds  Shear  Stess  Scaled  by  Local  Velocity  at  Height  from  Wall 
Y  =  0.66<5o 

Turbulent  kinetic  energy  (TKE)  was  calculated  as  described  in  Chapter  4,  using  the  rela¬ 
tionship  v'  k.  w '  to  compensate  for  the  inability  to  measure  w1  with  the  present  test  setup.  In 
Figure  5.44,  the  TKE  for  each  wedge  flow  rises  upon  nearing  the  interaction  and  quickly  returns 
to  near  upstream  values  in  roughly  1.550  at  the  measurement  height  Y  =  0.66<5o.  The  value  at 
the  interaction  for  both  the  5"  and  7°  wedge  flows  are  nearly  equal  but  the  10°  wedge  flow  has  a 
significant  rise  at  the  interaction.  This,  again,  supports  the  conclusion  that  the  10°  wedge  flow  is 
fundamentally  different  and  may  be  separated.  However,  there  is  no  LDV  evidence  to  support  the 
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conclusion  that  the  7°  wedge  flow  is  separated.  Note  that  the  order  of  increasing  TKE  upstream 
of  the  interaction  is  reversed  downstream  of  the  interaction  with  the  5°  wedge  flow  retaining  the 
highest  TKE  after  the  interaction. 
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Figure  5.44  Turbulent  Kinetic  Energy  at  Height  from  Wall  Y  =  O.66<50 


5.9.1  Fluctuating  Mach  Number.  As  stated  in  the  discussion  on  the  flat  plate  results,  the 
fluctuating  Mach  number  is  a  measure  of  the  validity  of  Morkovin’s  hypothesis.  The  fluctuating 
Mach  number  was  calculated  from  the  data  and  examined  to  ensure  that  the  value  never  approached 
M'  =  1.0.  In  all  cases  the  maximum  values  of  M'  never  exceeded  0.40.  Caution,  however,  should 
be  exercised  using  these  data  near  the  incoming  and  reflecting  shock  waves. 

5.9.2  Wall  Shear  Stress.  By  integrating  the  boundary  layer  momentum  equation  from 
the  wall  to  a  point  in  the  boundary  layer  a  distance  y  from  the  wall  (where  y  must  be  within  in  the 
boundary  layer),  it  can  be  shown  that  the  boundary  layer  shear  stress  should  tend  towards  the  wall 
shear  stress  (11).  With  the  assumption  that  pu =  0  and  no  blowing  or  suction,  the  boundary 
layer  momentum  equation  becomes 
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rp  dP 

Txy  -  +  -Q^y 


(5.1) 


Solving  for  rw  gives 


tt  - 

‘xy 


dP 


(5.2) 


This  is  called  the  Couette  Flow  assumption  and  is  valid  at  values  of  y+  <  500.  Couette 
flow  also  assumes  the  flow  is  unseparated.  The  pressure  gradients  used  for  this  study  in  Equation 
5.2  were  obtained  by  assuming  a  linear  distribution  of  pressure  between  the  static  pressure  mea¬ 
surements.  The  pressure  gradients  used  are  then  the  slopes  of  the  linear  segments  between  the 
static  pressure  measurements.  All  upstream  pressure  gradients  were  assumed  to  be  zero,  in  order 
to  avoid  differentiating  possible  sensor  noise  between  two  upstream  static  pressure  measurement 
points.  The  pressure  gradient  used  for  the  downstream  values  are  shown  in  Table  5.4. 


\0°Wedge 

7°  Wedge 

5  °W  edge 

— 1.32<50 

0 

0 

0 

1.32<50 

596121 

-11368 

-183461 

2.64<50 

-144284 

-298307 

-183461 

Table  5.4  Pressure  Gradients  in  Pa/m  Used  in  Calculation  of  Wall  Shear  Stress 


At  the  X  =  1.32<50  station  for  the  10°  wedge  flow,  the  Couette  flow  assumption  breaks  down, 
giving  a  nonsensical  value  of  tw.  Therefore,  that  data  point  has  been  discarded.  The  remaining 
wall  shear  stress  points  are  plotted  in  Figure  5.45.  This  figure  shows  a  700%  increase  in  the  wall 
shear  stress  over  the  upstream  values  (solid  and  dashed  lines  in  Figure  5.45)  for  the  5°  wedge  case. 
The  7°  wedge  flow  does  not  show  an  increase  in  the  wall  shear  stress  at  the  X  —  1.32<50  location 
because  of  the  small  pressure  gradient  measured  there.  However,  at  the  X  =  2.64<50  station,  the  wall 
shear  stress  is  1000%  larger  than  the  upstream  value.  Along  the  same  lines,  the  10°  wedge  flowfield 
experienced  a  670%  increase  in  wall  shear  stress.  The  wall  shear  stress  dramatically  increases  due 
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to  the  shock  wave  interaction  and  has  not  returned  to  the  upstream  value  at  the  last  measurement 


station.  The  measured  y+  for  the  lowest  measurement  from  the  wall  for  this  study  was  315  ±  5.0. 
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Figure  5.45  Comparisons  of  Wall  Shear  Stress 

The  corresponding  skin  friction  coefficient  chart  is  shown  in  Figure  5.46.  The  skin  friction 
coefficient  was  calculated  using  Equation  5.3  and  the  values  of  rw  shown  in  Figure  5.45.  p0  and  U0 
are  the  reference  quantities  defined  in  Table  5.2.  The  skin  friction  coefficient  follows  similar  trends 
as  the  wall  shear  stress,  showing  a  dramatic  increase  due  to  the  shock  wave  interaction. 
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Figure  5.46  Skin  Friction  Coefficients  for  the  Shock  Wave/Boundary  Layer  Interaction 


VI.  Conclusions  and  Recommendations 


6.1  Summary  and  Conclusions 

This  study  had  two  goals,  to  increase  the  understanding  of  a  shock  wave/boundary  layer  inter¬ 
action  on  a  flat  plate  and  to  provide  information  to  increase  the  turbulence  database  for  compressible 
flows.  To  this  end,  static  surface  pressures  were  measured  through  the  shock  wave/boundary  layer 
interaction  and  Laser  Doppler  Velocimetry  (LDV)  was  used  to  measure  the  mean  and  fluctuating 
flow  characteristics  of  the  shock  wave/boundary  layer  interaction.  The  flow  properties  derived 
from  these  measurements  included  velocity  and  Mach  number  profiles,  the  velocity  fluctuations, 
Reynolds  shear  stress,  turbulent  kinetic  energy,  wall  shear  stress,  and  the  skin  friction  coefficient. 

The  results  from  Chapter  5  allow  some  general  conclusions  to  be  drawn.  The  flat  plate  results 
from  the  present  study  match  well  with  the  data  from  Luker  (13)  and  Hale  (9)  and,  also,  with  Elena 
and  LaCharme  (7).  Additionally,  the  present  shock  wave/boundary  layer  interaction  data  match 
with  Rose  and  Johnson  (16)  with  the  exception  of  the  cross-stream  velocity  fluctuation.  The  cross¬ 
stream  velocity  measurements,  throughout  the  present  experiment,  may  not  be  valid  over  ta  1.0. 
This  is  caused  by  shock  waves  passing  through  the  traverse  profiles  leading  to  a  consistently  lower 
number  of  coincident  samples  farther  from  the  wall.  Examination  of  the  downstream  velocity 
profiles  reveals  that  velocity  does  not  change  appreciably  in  the  streamwise  direction  for  a  particular 
shock  strength.  However,  the  wall  shear  stress  and  skin  friction  coefficient  increase  over  500%  across 
the  interaction. 

The  shock  wave/boundary  layer  interaction  causes  significant  changes  in  the  downstream 
flow  structure.  The  velocity  fluctuations  have  increased,  particularly  near  the  wall.  The  stream- 
wise  velocity  fluctuations  increased,  for  example,  29%  for  the  5°  wedge  flow  near  the  wall.  The 
one  exception  to  this  was  the  10°  wedge  flow  which  showed  a  lower  streamwise  velocity  fluctua¬ 
tion  downstream  of  the  interaction.  It  is  believed  that  the  favorable  pressure  gradient  off  of  the 
interaction-induced  separation  region  dampens  out  these  fluctuations.  The  cross-stream  velocity 
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fluctuations  also  increased  across  the  interaction.  The  5°  wedge  flow  showed  an  approximate  37% 
increase  near  the  wall  when  comparing  upstream  and  downstream  measurements.  The  Reynolds 
shear  stress  increased  by  60%  near  the  wall  for  the  5°  wedge  flow.  The  turbulent  kinetic  energy 
increased  by  as  much  as  29%  across  the  interaction  for  the  5°  wedge  flow. 

The  strength  of  the  impinging  shock  also  has  significant  effects  on  the  properties  downstream 
of  the  interaction.  The  mean  velocity  and  Mach  number  defects  caused  by  the  impinging  shock  wave 
increased  with  shock  strength,  as  expected  from  inviscid  theory.  Most  of  the  turbulence  properties 
also  increased  with  shock  strength  near  the  interaction  with  the  exception  of  the  cross-stream 
velocity  fluctuation.  The  streamwise  velocity  fluctuations  increased  with  increasing  shock  strength 
but  recovered  to  upstream  values  quickly  with  the  exception  of  the  10°  wedge  flow  as  discussed 
above.  In  addition,  the  Reynolds  shear  stress  and  the  turbulent  kinetic  energy  increased  with 
increasing  shock  strength  near  the  interaction.  However,  the  cross-stream  velocity  fluctuations  do 
not  appear  to  depend  on  shock  strength  for  the  range  of  shock  strength  considered.  One  exception 
to  this,  however,  was  the  10°  wedge  flow  very  near  the  interaction,  which  showed  a  significant 
increase. 

The  lower  number  of  collected  coincident  bursts  increased  the  error  bounds  of  this  experiment 
over  the  studies  of  Luker  (13)  and  Hale  (9).  However,  the  error  increase  was  small  and  the  study 
is  still  considered  acceptable  with  the  error  bounds  stated  in  Appendix  B. 

6.2  Recommendations 

These  data  should  be  compared  to  CFD  solutions  using  established  turbulence  models,  similar 
to  the  study  done  by  Fick  (8).  However,  several  aspects  of  this  experiment  could  be  improved  to 
make  the  data  more  useful. 
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The  first  recommendation  is  to  use  a  different  shock  generating  scheme,  instead  of  the  wedge 
shock  generator.  The  wedge,  roughly  the  same  height  as  the  incoming  boundary  layer,  caused  an 
unexpectedly  large  flow  turning  angle,  complicating  comparison  to  numerical  analysis. 

To  increase  the  turbulence  data  base,  additional  streamwise  profiles  should  be  taken  closer 
to  the  wall  to  confirm  the  effects  near  the  wall  of  the  shock  wave/boundary  layer  interaction. 
Streamwise  profiles  are  an  excellent  method  of  determining  the  effect  of  both  shock  strength  and 
location  in  a  shock  wave/boundary  layer  interaction. 

The  olive  oil  seeder  mechanism  should  be  redesigned  to  give  a  more  even  distribution  of  seed 
through  the  boundary  layer.  It  appeared  during  the  experiment  that  the  seeder  holes  as  presently 
designed  created  dead  zones  in  the  boundary  layer  which  had  very  little  seed.  This  caused  the 
number  of  coincident  bursts  collected  to  drop.  A  seeder  with  a  continuous  vertical  slit  may  solve 
this  problem. 

The  final  recommendation  is  to  take  accurate  wall  shear  stress  measurements  over  the  shock 
wave/boundary  layer  interaction.  The  wall  shear  stress  is  the  determining  factor  in  the  magnitude 
of  the  friction  drag  and  heat  transfer  to  the  surface.  The  current  LDV  setup  was  limited  in  how 
close  it  could  get  to  the  wall,  making  it  difficult,  if  not  impossible,  to  accurately  determine  the 
wall  shear  stress.  The  Couette  flow  assumption  gives  an  idea  of  the  magnitude  but  this  assumption 
breaks  down  as  the  shock  strength  is  increased.  Obtaining  accurate  wall  shear  stress  measurements 
would  aid  in  the  analysis  of  the  data  presented  here  and  would  also  be  useful  for  further  research 
performed  in  this  facility. 
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Appendix  A.  Rose  and  Johnson  1975  Data 


A.l  Rose  and  Johnson  1975  Tabular  Data 


Table  A.l  Mean  Velocity  Profiles 


X  =  13.7  cm 

X  =  21.3  cm 

X  =  23.8  cm 

y  (cm) 

U(m/s) 

y(cm) 

U  (m/s) 

y(cm) 

U(m/s) 

0.22 
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0.27 

413 

0.13 
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0.34 

500 

0.40 
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0.26 

406 

0.52 
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0.53 

456 

0.37 

430 

0.74 

0.65 

475 

0.52 

449 

0.99 

556 

0.77 
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0.65 
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1.26 
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mm 
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1.50 

■  ' 

m 
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1.00 
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1.76 

■Eg 
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1.27 
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2.02 
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1.51 
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1.76 
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2.02 

559 

2.27 
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Table  A. 2  Streamwise  Velocity  Fluctuations 


X  =  13.7  cm 

X  =  21.3  cm 

X  =  23.8  cm 

y  (cm) 

y  (cm) 

E£1 

y(cm) 

Ea 

0.10 

0.0967 

0.26 

0.1390 

0.13 

0.1510 

0.15 

0.0900 

0.38 

0.1238 

0.26 

0.1310 

0.21 

0.0857 

0.51 

0.1024 

0.38 

0.1210 

0.25 

0.0752 

0.64 

0.0943 

0.51 

0.1130 

0.50 

0.0667 

0.76 

0.0824 

0.64 

0.0995 

0.76 

0.0567 

1.00 

0.0671 

0.76 

0.0857 

1.00 

0.0479 

1.25 

0.0498 

1.00 

0.0695 

1.26 

0.0433 

1.51 

0.0248 

1.25 

0.0530 

1.51 

0.0357 

1.51 

0.0288 

1.76 

0.0300 

1.77 

0.0172 

2.02 

0.0158 

2.02 

0.0126 

2.28 

0.0133 

2.28 

0.0074 

2.54 

0.0119 

2.54 

0.0083 
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Table  A. 3  Crosswise  Velocity  Fluctuations 


X  =  13.7  cm 

X  =  21.3  cm 

X  =  23.8  cm 

y  (cm) 

<«'> 

u 

y  (cm) 

.<y> 

U 

y  (cm) 

<y> 

V 

0.18 

0.0529 

0.26 

0.0348 

0.26 

0.0514 

0.26 

0.0533 

0.38 

0.0333 

0.38 

0.0414 

0.38 

0.0448 

0.51 

0.0381 

0.51 

0.0324 

0.51 

0.0414 

0.64 

0.0333 

0.64 

0.0376 

0.64 

0.0376 

0.76 

0.0271 

0.76 

0.0395 

0.76 

0.0343 

1.02 

0.0167 

1.02 

0.0238 

1.02 

0.0295 

1.27 

0.0162 

1.27 

0.0257 

1.27 

0.0229 

1.52 

0.0214 

1.52 

0.0124 

1.52 

0.0229 

1.77 

0.0152 

1.77 

0.0181 

2.03 

0.0152 

2.03 

0.0238 

2.29 

0.0181 

Table  A.4  Reynold’s  Normal  Stress 


X  = 

=  13.7  cm 

X  =  21.3  cm 

X  =  23.8  cm 

y  (cm) 

pU'2(N/7Tl2) 

y  (cm)  pU,2(N/m2) 

y  (cm) 

pU,2(N/m2) 

| MM 

■ 

2637 

0.142 

i&K 

IES9B 

2693 

0.204 

810 

0.500 

2376 

0.383 

2743 

0.259 

676 

0.628 

2168 

0.500 

2713 

0.500 

648 

0.756 

1853 

0.628 

2366 

0.756 

543 

1.00 

1608 

0.756 

2000 

1.000 

438 

1.27 

1000 

1.00 

1582 

1.27 

390 

1.53 

276 

1.27 

1087 

1.53 

371 

1.53 

352 

1.77 

190 

1.77 

133 

2.03 

57 

2.28 

38 
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Table  A. 5  Reynold’s  Shear  Stress 


X  = 

=  13.7  cm 

X  = 

=  23.8  cm 

y  (cm) 

pu'  v'  (N  /  m2) 

y  (cm) 

pu'v'(N/m 2) 

0.174 

120 

0.253 

316 

0.253 

125 

0.383 

327 

0.383 

139 

0.500 

414 

0.500 

155 

0.628 

477 

0.628 

157 

0.756 

420 

0.756 

145 

1.00 

423 

1.00 

116 

1.27 

220 

1.27 

84 

1.53 

91 

1.53 

73 

1.77 

31 

1.77 

32 

2.03 

13 

2.03 

25 

2.28 
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A. 2  Rose  and  Johnson  1975  Graphical  Data 


Figure  A.l  Mean  Velocity  Profiles 


Figure  A. 2  Streamwise  Velocity  Fluctuations 
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Figure  A. 3  Crosswise  Velocity  Fluctuations 
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Figure  A. 4  Reynolds  Normal  Stress 
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Appendix  B.  Error  Analysis 


The  error  analysis  for  this  experiment  follows  Luker’s  Appendix  B  ’’Error  Analysis”  (13) 
closely,  as  that  reference  calculated  all  the  required  error  statistics  for  the  present  test  setup  and 
wind  tunnel  combination.  The  only  error  values  that  have  changed  are  the  turbulence  statistics 
values.  These  values  have  changed  due  to  the  difference  in  the  number  of  coincident  bursts  collected. 
The  lowest  number  of  bursts  collected  during  the  experiment  was  121.  This  value  will  be  used  to 
determine  the  effect  on  the  error  in  the  turbulence  statistics,  to  highlight  the  worst  case  .  All  other 
variables  continue  to  have  the  same  error  values  as  presented  in  Luker. 

Define  the  error  from  the  low  number  of  bursts  as  eg  and  the  number  of  bursts  collected  as 
Nb  ■  Then  the  error  in  the  Root  Mean  Square  (RMS)  values  is  given  by 

«  =  100|«B-1-NB|  =  100  =  0  833%  (B.U 

A  b  Nb 

This  error  is  then  combined  with  the  errors  found  by  Luker  using  the  Euclidian  norm  of  the 
independent  error  terms.  For  example,  the  total  error  eu'Total  is  given  by 


(B.2) 


Using  equation  B.2,  £u'Tot<xl  —  1.89%.  The  remaining  error  terms,  both  Luker’s  error  cal¬ 
culations  and  the  present  study’s,  are  shown  in  table  B.l.  The  different  values  for  this  study  are 
marked  with  an  asterisk. 
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Appendix  C.  Turbulence  Statistics 


C.l  Skewness 


The  skewness  was  computed  using  the  method  described  in  Chapter  4  and  is  presented  below 
in  figures  C.l  through  C.16.  The  error  bounds  for  these  data  are  presented  in  appendix  B  and  the 
data  is  tabulated  in  appendix  D. 


2.5 

1  1  1  i 

.  | 

2.0 

■ 

: 

■  A  ♦ 

■  ♦ 

A 

o 

■ 

A  ♦ 

i©°  1.5 

- 

O  *  ■ 

- 

O 

♦  ■  A 

O 

♦  A  ■ 

1.0 

- 

- 

■  AO 

O  ■  A 

■  5°  Wedge 

- 

A  7°  Wedge 

«C 6 

0.5 

- 

♦  10°  Wedge 

- 

- 

OFIat  Plate 

■ 

.  .  .  1 

,  1  ,  .  ,  ,  ° 

-15 

-10 

-5  0 

U  Skewness 

Figure  C.l  X  =  —1.32 60  Streamwise  Skewness 


C.2  Flatness 

The  flatness  was  computed  using  the  method  described  in  Chapter  4  and  is  presented  below 
in  figures  C.17  through  C.32.  The  error  bounds  for  these  data  are  presented  in  appendix  B  and 
the  data  is  tabulated  in  appendix  D. 
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Figure  C.3  X  =  Q60  Streamwise  Skewness 
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Figure  C.5  X  —  1.3250  Streamwise  Skewness 


Figure  C.6  X  =  1.32<$0  Cross-stream  Skewness 
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Figure  C.ll  5°  Wedge  Streamwise  Skewness  at  each  Station 
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Figure  C.15  10°  Wedge  Streamwise  Skewness  at  each  Station 


Figure  C.17  X  =  — 1.32<50  Streamwise  Flatness 


Figure  C.26  Cross-stream  Flatness  at  a  Height  from  Wall  Y  =  0.665, 


Figure  C.27  5°  Wedge  Streamwise  Flatness  at  each  Station 


Figure  C.32  10°  Wedge  Cross-stream  Flatness  at  each  Station 


C.3  Intermittency 

The  intermittency  was  computed  from  the  flatness  using  the  method  described  in  Chapter  4 
and  is  presented  below  in  figures  C.33  through  C.39.  The  error  bounds  for  these  data  are  presented 
in  appendix  B  and  the  data  is  tabulated  in  appendix  D. 


C-4  Velocity  Correlation  Coefficient 

The  velocity  correlation  coefficient  was  computed  using  the  method  described  in  Chapter  4 
and  is  presented  below  in  figures  C.40  through  C.46.  The  error  bounds  for  these  data  are  presented 
in  appendix  B  and  the  data  is  tabulated  in  appendix  D. 
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Figure  C.33  X  =  —1.32 60  Intermittency  Function 
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Figure  C.36 
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Figure  C.37  5°  Wedge  Intermittency  Function 
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Figure  C.39  10°  Wedge  Intermittency  Function 


Figure  C.41  Velocity  Correlation  Coefficient  at  X  =  05, 


Figure  C.42  Velocity  Correlation  Coefficient  at  X 


Figure  C.44  5°  Wedge  Correlation  Coefficient 


Appendix  D.  Tabulated  Data 


The  tabulated  data  for  all  coincident  data  points  is  presented  in  Figures  D.l  through  D.24. 
The  density,  p,  has  units  of  ^  while  the  velocities  have  units  of  ™ .  The  temperature  is  in  degrees 
Kelvin. 


D.l  Flat  Plate  Tabulated  Data 
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Figure  D.l  Tabulated  Data  for  Flat  Plate  at  x  =  48  cm 
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D.2  5°  Wedge  (9eff  =  8.2°)  Tabulated  Data 


Figure  D.2  5°  Wedge  ( 9ejj  =  8.2°)  Flow  X  =  —1.32 50  Coincident  Data 
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Figure  D.5  5°  Wedge  (Oefj  —  S.2°)  Flow  X  =  2.64 60  Coincident  Data 
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Figure  D.6  5°  Wedge  (Oeff  —  8.2°)  Flow  X  =  2.64<50  Coincident  Data  (Cont.) 
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Figure  D.7  5°  Wedge  (Oejj  =  8.2°)  Flow  Y  =  0.66So  Coincident  Data 
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Figure  D.8  5°  Wedge  ( 9eff  =  8.2°)  Flow  Y  =  0.666o  Coincident  Data  (Cont.) 


D.3  7°  Wedge  (6efj  =  10.4°^  Flow  Tabulated  Data 
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Figure  D.9  7°  Wedge  (0eff  =  10.4°)  Flow  X  =  — 1.32<50  Coincident  Data 
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Figure  D.10  7°  Wedge  (Oeff  =  10.4°)  Flow  X  =  0 60  Coincident  Data 
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Figure  D.ll  7°  Wedge  (0e//  =  10.4°)  Flow  X  =  1 .32(50  Coincident  Data 
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Figure  D.12  7°  Wedge  ( 9ejj  —  10.4°)  Flow  X  =  1.32 5a  Coincident  Data  (Cont.) 
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Figure  D.13  7°  Wedge  (9ejj  =  10.4°)  Flow  X  =  2.64<50  Coincident  Data 
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Figure  D.14  7°  Wedge  (9e/f  =  10.4°)  Flow  X  =  2.64<50  Coincident  Data  (Cont.) 
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Figure  D.15  7°  Wedge  (0e//  =  10.4°)  Flow  Y  =  0.66<5„  Coincident  Data 


D-15 


X/5 

v/uc 

vyu 

VAJC 

SKV 

FLy 

Bifmasn 

KIM  MUS  MO* 

Miin«i:iain« 

BIMI*W*>« 

B.UUIIJiM 

BiiMiioimia 

EiUM&l 

EESS33S3 

BiVW.<VJ:l>iB 

B!HK'iTT!T!TI 

tmmm 

mtvmzm 

w&mite urn 

■Hdftf.MglDK 

BiU'.’AMJiW 

BataoistM 

0.7692  IE-03 

0.20584E-01 

0.I9949E-01 

B'H&lhiM 

-0.75000E+00 

0.2343  8E-02 

0.21735E-01 

0.14161E+02 

0.384 54E-03 

Bairmarsa 

0.33775E-03 

K35SSS! 

BiHijlAIJilB 

BUEaa 

O.OOOOOE+OO 

-0.45689E-01 

BiWAIMigilB 

0.98131E+01 

B4Hili»UiVB 

0.25000E+00 

-0.S5263E-01 

IMl'MM 

0.49320E+00 

an 

0.50000E+00 

-0.4974  IE-01 

0.26764E-0I 

BIMMiMiW 

0.72639E+01 

0.75000E+00 

-0.47962E-01 

»!?»:»=», IB 

0.11415E+02 

-0.49274E-01 

BiMIMIJiW 

EEEBSa 

■ly^waiTB 

-0.42594E-01 

0.26849E-01 

Bil>jkL'!l  JilB 

■iUiIMJMB 

0.15000E+01 

-0.38568E-01 

0.24920E-01 

•  ■ 

■illEKIgfiTM 

CE2EHS' 

0.17500E+0! 

-0.31207E-01 

0.24007E-01 

■innh^iT 

0.200 00E +01 

-0.24438E-01 

MiWWAJtnw 

■igEiUtiaflil 

Miutian 

■Hwinmrcnf 

-0.19435E-01 

MiW.WIigilB 

BiMHrflAM 

■lEiriMIJilM 

HSE353I 

BEHESl’iB 

-iMl'.'i'.iJiM 

OEUESEJ 

liHtliMM 

1  0.27500E+01 

IB+HMJiM 

iTB 

wMm&m 

■iWiMHJi.18: 

-0.81408E-02 

0.24047E-0! 

W 

■iHfcttlllsRlIlB 

<*B®i 

X/8 

Ruv 

TKE 

y+ 

lEMmi 

ggrogiai 

BilWMi&Mi 

-0.17500E+01 

EEB^3I 

■0EElK)iffi!iB 

BHIE5353 

WW8WBW 

■HkT.r.  :-i7B 

BilkijirUiVB 

Bnwn:i):riHB 

■fr.Tr^iaty 

■m  wif.maiEB 

bsessem 

eaiiMiMfliii 

MM.’MWrm 

Ksa 

E33SS&H 

mikkmmim 

Mmmxm 

Kusmmm 

BB5553S51 

Baaaiaai 

KiygMO-riTil 

EEZMga 

BiUTBtlJiM 

■niy.mniMnif 

KUEWtiMMtfll 

mum 

KEjUBga 

kixi:i:m:i]gfliB 

■uraiiiiiiafTil 

ihseisjS 

BB2335M 

Kmcir.iUMtEB 

■uw/yjxgUM 

■*<n«i«i*B 

B5EES35fli 

■niHMsa 

■»w«!Maa!il 

BilUKW 

assaaga 

KESSSSM 

BE 

■snasgs 

EE^35E 

XiMUdSMW 

■MZiwginil 

BES5B5M 

BilliWrfaJil»E 

bih«ibh»SK 

KH1E3353I 

Hm 

;  EEgffi 

■tiwagRiaan'l 

BjSEEiigffil 

EEBBS3I 

E  Ha 

nmisHMCia 

:  '  MMF..- 

KBS5353 

RKn 

BiHiiikVi 

■iHd-HHaW.l 

■itBMiiB 

MiJIWiWXT&W 

rnamdsm 

■nikviaTTB 

wmamm 

BilllltfUiW 

BBSi 

■lIlk.'.-JlIB 

mnmjrim 

Figure  D.16  7°  Wedge  (0e//  =  10.4°)  Flow  Y  —  0.666„  Coincident  Data  (Cont.) 


D-4  10°  Wedge  (9ejj  =  13. 5°)  Flow  Tabulated  Data 
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Figure  D.17  10°  Wedge  ( 9ejj  =  13.5°)  Flow  X  =  — 1.32<50  Coincident  Data 
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Figure  D.18  10°  Wedge  ( 0ejf  =  13.5°)  Flow  X  =  0So  Coincident  Data 
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Figure  D.19  10°  Wedge  (8eff 


13.5°)  Flow  X  =  1.325,,  Coincident  Data 


D-19 


Y 

V/Ue 

wu 

V'/Ue 

SKV 

FLV 

tjy/pu2 

■am 

0.58523E-01 

EESBUEM 

BUS! 

0.30000E+01 

wmuwm 

KH 

Ej^33B5B 

EEB3BSM 

0.40000E+0 1 

BBaaagai 

laaaBaaa 

■iMwaW 

■iTKBSEBM 

■imiuhhuciiT 

-0.84526E-01 

■iMAvtigilB 

BBEH3EM 

jWBSMSai 

Eggai 

eikcimmiuib 

0.90000E+01 

EMI 

0.10000E+02 

WSnEBEi 

0.257 81  E-01 

0.22127E+02 

0.1  IOOOE+02 

-0.12195E+00 

Tsmmm 

EHS 

0.12000E+02 

-0.1 1788E+00 

■  m 

0.49668E+02 

-0.1 1406E+00 

gffw.oansngil 

EM3Sa 

0.14000E+02 

-0.10948E+00 

WMEBU 

■lEUPmgriTl 

lilrf.nltfcVrrfPM 

0.15000E+02 

-0.10703E+00 

0.28I51E-01 

wmmm m 

-0.26869E+01 

0.16000E+02 

-0.10225E+00 

I.MBM 

■OMr.maiU 

EE3ESI 

0.17000E+02 

-0 .97044 E-01 

■iHcMdltgilB 

-0.46237E+01 

O.4646OE+02 

■lU’UIJiW 

0.18000E+02 

-0.926 15E-01 

ESESE33B 

■inrisr.T.idfiTM 

EEsasa 

EHE5353 

0.19000E+02 

ME3EE:® 

gnE»^ish?ni 

K5EE5jgi3 

gjOTjTCibHiM 

gnnicii:.igiTM 

**»*■ 

BIEUMAIB 

BBJIHai!! 

Y 

> 

cc 

TKE 

y* 

0.20000E+01 

■iMBwaaifW 

ezshsh 

■iiiitKiaif 

0.31296E+03 

■mMtoiiKiTB 

im a 

BwaaitW 

B^aaa 

BHE5BS3 

WEBSEM 

■inittlilimitfiTW 

0.I8019E-02 

EME5SJ 

BHESHa 

EISSEEl 

■wjnwtisamM 

gattViuitilB 

0.20426E+00 

0.18783E-02 

■OWMIlgfiSl 

■i»r.<iiii:<i^nT 

■iwmiaiM 

fcie»s».iaiw 

Ezmsa 

BSESS3SJM 

■33ES2HB 

'  t 

wxLW&>m 

wunitur.iTm 

HiJmMigiaross 

BEM3SB 

wukrz m 

KMl 

SESEMiM 

■a’niiiimgiuM 

munn«m 

hiwm<rfll 

■aatgwaJ»BI 

BM-'lil.IrfMM 

KaasaaaaBi 

WSEEMEM 

wsmssma 

BBT  *  v 

■iiiiOTaif 

■i  KfttelXU'/M 

BHES3S3 

B  •  .  ' 

EBSESEl 

■ominmisrfiEB 

BESSIE 

EfltttiaEiI 

■  r‘  _ 

MiXiUMzF.im 

EESS1EI 

■tdKIMlaniEM 

E5S3ZUi3 

wnmzmm 

BESTHEJ 

wmmzsm 

■mumnistfiEB 

■3ESI2EH 

mnm#M 

MIM3M 

TFT 

EEESSEESi 

wni’Mm'rm 

■WfflfliM 

BESE2!I 

Wwnw.m 

1352333531 

-0.9438  IE-05 

■SHFHaiS* 

EBEB3M 

MiW.w-mm 

■iimtitinaroTB 

mmssssM 

KiilMJil 

BESEBiEB 

EBSEjSEl 

BMMIrJsRiEI 

■HEtHlia'Tl 

EB5EE1 

fcmmigiil! 

BEE5DEEB 

EE0H323I 

EB535J 

Biirai=!i!l 

■iHli'rtl:*.. 

■towikiimua 

msmmm 

■HKsMffli 

Figure  D.20  10°  Wedge  (0e//  —  13.5°)  Flow  X  =  1.32<50  Coincident  Data  (Cont.) 
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Figure  D.21  10°  Wedge  (6eff  =  13.5°)  Flow  X  =  2.64 60  Coincident  Data 
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Figure  D.22  10°  Wedge  (9eff  =  13.5°)  Flow  X  =  2.64(50  Coincident  Data  (Cont.) 
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Figure  D.24  10°  Wedge  (9ejf  =  13.5°)  Flow  Y  =  0.666o  Coincident  Data  (Cont.) 
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